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ABSTRACT 


The major differences in characteristics of propellers for air- 
craft gas turbines as they affect propeller design are outlined and 
discussed. The higher powers that appear probable with the 
aircraft gas turbine make for airplane speeds in the range in which 
it has been predicted that good propeller efficiency could not be 
obtained, but it is shown that good propeller efficiency can be 
maintained up into the range of 600 to700m.p.h. The configura- 
tion, of propellers for this type of application is discussed. An- 
other characteristic of the gas turbine is that it is essentially a 
constant-speed power plant which introduces new considerations 
in the design of propellers to give maximum efficiency at high 
speed and cruising and at take-off. Effect of diameter and 
- solidity on this problem are considered and the advantages of 
two-speed reduction gearing are discussed. This paper considers 
single rotation propellers primarily but dual rotation types are 
examined briefly. A comparison of thrust and thrust horsepower 
for turbojet and turbine-propeller power plants is made and the 
results are applied to a fighter-type airplane. Special control 
characteristics of turbine propellers and weight and vibration 
problems are also discussed. 


INTRODUCTION 


T° OBTAIN MAXIMUM PERFORMANCE from propellers 
for aircraft gas turbines, the propeller design must 
take into account several important characteristics of 
this type of power plant which differentiate it from the 
conventional reciprocating aircraft engine. These char- 
acteristics are: 

(1) Powers are higher, making higher airplane 
speeds possible. Whereas speeds of 450-500 m.p.h. 
have been the maximum obtained with present piston 
engines, the larger turbines should make speeds of the 
order of 600 m.p.h. and up possible. 

(2) The turbine is essentially a constant-speed 
power plant and gives lowest specific fuel consumption 
at maximum power. It will, therefore, tend to be oper- 
ated at much higher cruising r.p.m.’s and powers than 

the piston engine. 
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(3) Overspeed and overtemperature limits are more 
critical for the turbine and, hence, propeller control 
must be more precise. Present indications are that the 
propeller control and the turbine control should be 
combined for best results. 

This paper will discuss the effect of these three char- 
acteristics of the aircraft gas turbine on propeller design. 
Data are presented also to compare the performance of 
turbine-propeller and turbojet installations. 


AERODYNAMICS FOR HIGH SPEED 


Many predictions were made several years ago both 
within and without the propeller industry that catas- 
trophic reductions in propeller efficiency would be ex- 
perienced at speeds of about 400 m.p.h. Ata later date, 
when fighter airplanes had exceeded 400 m.p.h., the 
critical speed was revised to 500 m.p.h. To date, pro- 
peller-driven airplanes have attained level flight speeds 
of slightly over 500 m.p.h. with good propeller efficiency 
being maintained. In achieving this result it has been 
found that blade sections should be thin and that well- 
faired blade shanks are necessary, but that otherwise 
a conventional configuration is satisfactory. 

This experience indicates that propellers for gas tur- 
bines for airplanes in the 400-500 m.p.h. category need 
not be appreciably different from those for piston en- 
gines of the same rating so.far as aerodynamics are 
concerned, although later in the paper the effect of con- 
stant r.p.m. operation on the propeller characteristics 
is discussed. Fig. 1 shows the comparative perform- 
ance of a propeller for both gas-turbine and turbosuper- 
charged piston engine types of power plants for a large 
airplane having a maximum speed of 400 m.p.h. at 
30,000 ft. It will be noted that the same power rating 
at altitude is assumed, and that at sea level the take-off 
thrust for the turbine propeller is appreciably greater 
than for the piston-engine propeller even though the 
thrust in terms of pounds per horsepower is less. This 
figure also shows the advantageous effect of a higher 
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Variation of thrust with airplane velocity for piston 
engine and gas turbine. 


Fic. 1. 


propeller speed such as would be obtained with a two- 
speed reduction gear. 

For higher speeds of the order of 600 m.p.h. the prob- 
lem changes appreciably. In considering speeds of this 
magnitude it is desirable first to examine the location 
in the Mach Number spectrum of both the airplane and 
the propeller. For example, for the 600-m.p.h. airplane 
the airplane Mach Number is 0.785 at sea level and 0.88 
at 30,000 ft. With a rotational tip speed (xn D) of the 
propeller of 800 ft. per sec. the tip Mach Numbers for 
the two altitudes are 1.06 and 1.19, respectively, and 
the shank Mach Numbers for a propeller of conven- 
tional proportions having a spinner-to-propeller diame- 
ter ratio of between 0.15 and 0.2 are 0.80 and 0.90. 
Figs. 2 and 3 present these data for speeds from 500 to 
700 m.p.h. 

To determine the efficiency obtainable at these higher 
speeds, it is of interest to look at the effect on the ef- 
ficiency of the L/D ratio. Considering the propeller 
blade section at 0.75R as representative of the blade, 
the wind angle for 600 m.p.h. forward speed and 800 
ft. per sec. rotational tip speed is 55° 45’. The ef- 
ficiency of this section based on the blade element theory 
is 89.7 per cent for an L/D ratio of 20. Since this is a 
higher L/D ratio than can be expected in the-sonic re- 
gion, the effect of lower ratios is of interest. From avail- 
able wind-tunnel data, it is believed that L/D ratios of 
the order of 10 can be obtained throughout the trans- 
sonic range, which results in an efficiency for the 0.75R 
section of 79.5 per cent. Even if the ratio drops to 7 
the possible efficiency is 72.5. Hence, there is every 
reason to believe that we can expect reasonable ef- 
ficiency at speeds in the 600-700 m.p.h. range. 
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FLIGHT TESTS AT HIGH SPEED 


In considering what full-scale data are available to 
support the above conclusion, the tests conducted on the 
P-47] by Republic Aviation Corporation are of interest. 
The information released by the Army Air Forces 
shows that this airplane did slightly over 500 m.p.h. at 
about 34,000 ft. The propeller was a Curtiss C642S-B 
with 836-6C2-18 (13 ft. 0 in.) blades. Fig. 4 gives the 
data on this propeller and also the speed and Mach 
Number for each section of the blade. It will be seen 
that the tip Mach Number was 1.11 and that most of 
the blade was working appreciably above the so-called 
“critical’”” Mach Number. A strip analysis for this 
propeller for the flight-test condition gives a calculated 
efficiency of 74 per cent. Because of uncertainty as to 
the amount of exhaust jet thrust available from the 
turbosupercharger and the effect of compressibility at 
the test velocity on the airplane drag coefficient, it has 
not been possible to obtain an efficiency value that is 
considered reliable by back figuring from the known 
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Fic. 5. A comparison of experimental and calculated efficiency 
vs. airplane Mach Number. 


data, but from such checks as can be made it is believed 
that the strip analysis data are not far from correct. 

A further test of interest is that conducted by the 
Curtiss-Wright Propeller Division using an A.A.F.- 
supplied P-47D-30 airplane equipped with thrust and 
torque meters. By diving the airplane it was possible 
to reach extremely high forward Mach Numbers and 
to make thrust measurements from which shaft thrust 
efficiency could be calculated. Fig. 5 presents the re- 
sults of one of the tests that was conducted using a 
standard production propeller with 836-14C2-18 blades. 
An efficiency of over 80 per cent is maintained up to 
500 m.p.h. At the maximum speed of 592 m.p.h. 
reached in this test, the airplane Mach Number being 
0.84, the efficiency was 62 per cent. These values are 
considered particularly significant because they are for 
a propeller designed for 450 m.p.h., operating at con- 
siderably less than the design lift coefficient of its blade 
sections. 

Since the above results were obtained with present 
conventional propellers designed for speeds of about 450 
m.p.h. and, since there are a number of improvements 
in design which, as discussed later in the paper, appear 
certain to increase the efficiency to at least 70 per cent, 
it is evident that full-scale data support the conclusion 
that good efficiencies can be obtained at speeds in the 
600 m.p.h. and over range. 


DESIGN CHANGES TO INCREASE EFFICIENCY 


In considering what can be accomplished by changes 
in propeller design to increase efficiency in high-speed 
flight, two possibilities are of particular importance. 
First, is there any development, either conventional or 
unconventional, which will increase the L/D ratio in 
the sonic speed region? If there is, the efficiency will 
be increased. Second, can we design a propeller in 
which only that part of the blade that produces maxi- 
mum thrust is exposed? If this is possible, the problem 
of shank drag (which, when large, reduces propeller 
efficiency) will disappear. 

Thinner sections are one development that will aid 
in increasing blade section L/D and in reducing shank 
drag. The trend for some time has been toward the 
use of thinner sections, and it has been found practical 
with the newer types of blade construction and the 
better knowledge of vibration conditions to make sub- 
stantial thickness reductions. Before World War II, 
thickness ratios were usually 0.075 or greater at the 0.75 
radius but they had been reduced to 0.055-0.065 by the 
end of the war. Shank sections have also been sub- 
stantially reduced in both thickness ratio and in abso- 
lute thickness during this period. 

More efficient airfoil sections would also offer further 
possibilities in the direction of improved efficiency but 
there is no information to indicate that any section bet- 
ter than the NACA 16 series will be available in the near 
future. The possibilities in this direction are far from 
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exhausted, however, and more research work is strongly 
indicated. It might be well to mention at this point 


that despite all the wind-tunnel work that has been 


AERONAUTICAL 








SCIENCES—SEPTEMBER, 1946 





TABLE 4 
Data on Propeller Blades for 1,000 HP. at Various Speeds 











V—M.P.H. 500 600 700 
Hp. 1,000 1,000 1,000 1,000 1,000 
«nD 800 800 800 900 900 
Diameter 13 ft. Oin. 13 ft.O0in. 13 ft.Oin. 13 ft. Oin. 11 ft. Oin 
A.F. 135 100 75 63 90 
Tip Mach Num- 
ber 1.065 1.169 1.279 1.34 1.34 


Maximum design 


CL .55 0.55 0.55 0.55 0.55 
h/b—0.75R 0.05-0.06 0.05-0.06 0.05-0.06 0.05-0.06 0.05-0. 0¢ 
h/b—spinner 

junction 0.12 0.12 0.12 0.12 0.12 

81.5 78.5 72.5 73.1 74.5 


»—— Per cent 





done, there is still a dearth of accurate data on any air 
foils in the transsonic range. 


Unconventional Configurations 


Another possible way of increasing efficiency at high 
speeds is to apply the principle of sweepback to a pro- 
Fig. 6 shows a sketch of a conventional 
At 600 m.p.h. the efficiencies 


peller blade. 
and a swept-back blade. 
calculated by strip analysis show 75.5 per cent for the 
conventional plan form and 82.0 per cent for the swept- 
The operating conditions used in the 
Incidentally, this is 


back type. 
analysis are given on the figure. 
the only important development on propellers which 
has come from Germany. Fig. 7 is a replot of data 
given in a paper entitled ‘“‘Aerodynamische und Flug- 
mechanische Fragen der Luftschraubenentwicklung,” 
by August Wilhelm Quick, which was published in the 
Schriften der Deutschen Akademie der Luftfahrtforschung 
for April 2, 1943. This curve is believed to be based 
on propeller model tests and shows an improvement 
with sweepback of about 7 per cent at the highest tip 
Mach Number. No full-scale data are yet available 
on swept-back blades but it is hoped to have flight-test 
results within the next few months. 

The unconventional configuration of a blade in which 
only the thrust-producing part of the blade is exposed 
and all of the drag-producing part of the shank is housed 
presents some interesting possibilities. If we assume 
a spinner diameter of 3 ft., then the calculated efficiency 
for a single blade absorbing 1,000 h.p. is 81.5 per cent 
for 500 m.p.h., 78.5 per cent for 600 m.p.h., and 74.5 
per cent for 700 m.p.h. The data for these propellers 
are given in Table 1. The practical application for 
such a configuration will probably be as a relatively 
small-diameter propeller using eight or ten blades ap- 
proximately as shown in Fig. 8. Diameters of 10-12 
ft. for 5,000 h.p. appear possible with this type of pro 
peller. 

A further promising type of propeller is one in which 
the entire exposed portion of the blades operates at 
supersonic speed. The object of this arrangement is 
to have all of the working part of the blade operating 
above about 1.4 Mach Number in order to get above 
the speed at which the high drag coefficients associated 
with the transsonic region are experienced. A further 
reason is to keep the blade angle down to about 45° 
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where maximum efficiency at any given L/D is obtained. 
A propeller of this type would be similar to that shown 
in Fig. 8 except that the diameter would be smaller. 

A somewhat secondary but still important factor in 
designing a propeller for a turbine is to provide for 
maximum pressure recovery for the turbine intake air. 
This, in general, means a minimum blade shank thick- 
ness ratio which can best be achieved by the smallest 
thickness per blade and the minimum number of blades. 
This would indicate that for conventional types of pro- 
pellers three blades would be preferable if adequate 
total solidity can be obtained. The fairing of blade 
shank to spinner is important and must permit pitch 
change and, at the same time, provide an aerodynami- 
cally clean junction at all blade angles in the operating 
range. 

Another configuration that gives promise is that 
known variously as the ‘“‘annular’’ spinner, blower 
spinner, or NACA Type E cowl. This arrangement is 
shown on Fig. 8 as applied te the multiblade type of 
propeller, but is equally applicable to a conventional 
propeller of three or four blades. Because the air is 
slightly slowed down before it passes the blade shanks, 
there appears to be a good possibility of greater pres- 
sure recovery and, hence, more turbine power. 


CRUISING AND TAKE-OFF PERFORMANCE 


As stated in the introduction, the gas turbine is most 
efficient at maximum power where it has its lowest 
specific fuel consumption. Hence, it will tend to be 
operated for cruising at close to maximum rating. In 
order to determine the probable size of propellers that 
will be used for turbine power plants, an analysis has 
been made for a turbine having a cruising rating of 3,000 
h.p. at 25,000 ft. and a take-off rating of 5,400 h.p. at 
sea level. Fig. 9 gives the efficiency for cruising at 400 
m.p.h. for propellers of diameters of 16 to 22 ft. but of 
the same solidity. Fig. 10 shows the same data for 
propellers of 18-ft. diameter and several solidities. 
Fig. 11 shows the static thrust for these same propellers. 

These data are presented in Table 2 where it will be 
noted that if the r.p.m. is chosen for peak cruising ef- 
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Fic. 10. Propulsive efficiency vs. propeller r.p.m. for cruise at 
400 m.p.h. at 3,000 hp. at 25,000 ft. 
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Static thrust vs. propeller r.p.m. for T.O. rating of 
5,450-hp. for three-blade propellers. 


DIAMETER ACTIVITY FACTOR 
160" 120 
1s'o" 80 
1s'o 120 
18’o" 160 
20'0 120 
22'0" 120 





“TABLE 2 


Propulsive Efficiency and Static Thrust for Propellers of Varying Diameter and Solidity at Various R.P.M.’s 
Cruising Rating—3,000 Hp., at 400 M.P.H. at 25,000 Ft. 


Take-Off Rating—5,400 Hp. at 0 M.P.H. at Sea Level 





a4 tatceursdhlls ‘Dual 
—_—_—_—_—_—__——— Single Rotation—_—___-__—___. Rotation 
Diameter 16 ft. 18 ft. 20 ft. 22 ft. 16 ft. 
A.F. 120 80 120 160 120 120 120 
No. Blades 3 3 3 3 3 3 6 
R.p.m. for maximum 7 n 80.7 81.0 82.2 81.9 82.6 82.7 83.0 
R.p.m. 890 795 785 760 690 610 820 
Ty 7,000 6,800 8,800 11,000 10,200 11,500 13,000 
R.p.m, for maximum 7 + 50 r.p.m. n 79.9 80.8 81.8 81.7 82.2 82.2 82.6 
R.p.m. 940 845 835 810 740 660 870 
To 7,900 7,600 9,700 12,300 11,500 13,200 13,500 
R.p.m. for maximum 7 and Ty, at n 80.0 81.0 82.2 81.9 82.6 82.7 83.0 
1.3 X this r.p.m, R.p.m. 1,160 1,030 1,020 990 900 795 1065 
To 12,100 10,600 13,800 17,200 16,100 18,200 16,200 
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ficiency for the various propellers the static thrusts are 
considerably below the maximum. [If slightly higher 
propeller r.p.m.’s are selected, then the take-off thrust 
is increased. The data in Table 2 show the effect of a 
50 r.p.m. increase and of a take-off r.p.m. 30 per cent 
higher than that for peak efficiency as could be obtained 
with a two-speed reduction gear. 

A similar analysis for a six-blade dual rotation pro- 
peller has been made and the results are also given in 
Table 2. The data presented in Figs. 9-12 and in 
Table 2 can be used for other powers by scaling up or 
down in the ratio of the square of the diameter. 
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Another approach to the problem of obtaining good 
cruising performance is the “double turbine” arrange- 
ment in which two turbines-drive one propeller. For 
long-range cruising, one turbine is declutched and the 
other is operated at its maximum continuous and, 
hence, most efficient rating. 


COMPARISON WITH TURBOJET 


It is interesting ‘to compare the relative performance 
of turbojet and turbine-propeller power plants. This 
can be done on several bases. One basis is the same 
rate of fuel flow for the two types, which is logical from 
two standpoints. First, the same rate of fuel flow rep- 
resents the same rate of potential power, that is, b.t.u.’s 
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per hour going into the power plant, and, hence, pro- 
vides a good measure of efficiency. Second, at any 
given stage in gas-turbine development it should be pos- 
sible to build a certain maximum size of compressor. 
This could then be used in either a turbojet or turbo- 
propeller unit and, again, a good measure of relative ef- 
ficiency can be obtained. 

The figures presented here are the results of a study 
made on this hypothesis of the same fuel consumption. 
Available data on the two types of units were scaled 
up or down in order to arrive at the basic figures of 
comparison. The propeller used was purposely chosen 
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Fic. 18. Turbine-propeller vs. turbojet (same fuel consumption). 
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Fic.19. Turbine-propeller vs. turbojet (same fuel consumption). 


25,000 ft. 


as one of small diameter and conservative assumptions 
as to the propeller efficiency above 500 m.p.h. were 
made. Figs. 12-17 show the thrust horsepower and 
thrust comparisons for altitudes of sea level, 25,000, 
and 45,000 ft. The superiority of the turbine-propeller 
is marked at all except the absolute maximum speed 
and, even then, the propeller is slightly superior. 

These data were then applied to a fighter-type air- 
plane. This type was selected because it is the class in 
which it is thought by many people that the turbojet 
will certainly supersede the propeller. To be conserva- 
tive, 2,000 Ibs. increase in gross weight and 5 per cent 
increase in drag were assumed for the turbine-propeller 
airplane. Figs. 18 and 20 show that the overall per- 
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Fic. 20. Turbine-propeller vs. turbojet (same fuel consumption). 
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formance of the propeller-driven airplane is substan- 
tially superior to the turbojet. Maximum speed is 
about the same, but the propeller fighter, despite its 
higher weight, has a rate of climb nearly double that of 
the jet. Take-off run will be approximately one-half 
that of the jet-powered airplane. 

A calculation as to range shows that the propeller 
fighter has a definite advantage and for endurance at 
lower speeds the margin is large, the improvement for 
the propeller-driven fighter at 400 m.p.h. at 25,000 ft. 
being about 40 per cent. Another characteristic of 
great interest from the operational standpoint is that 
the propeller fighter at 350 m.p.h. at 25,000 ft. has a 
turning radius of 3,200 ft. as compared with 4,400 ft. 
for the turbojet fighter. 

It might be argued that the probable future improve- 
ments in fuel economy or momentary power increases 
of turbojets will make the present comparison invalid. 
This is not believed a proper objection, however, inas- 
much as any improvement in efficiency due to internal 
improvements in the turbojet will be equally appli- 
cable to the propeller turbine and, hence, the same rela- 
tive position will be maintained. 


WEIGHT AND VIBRATION 


Propeller weight for any given size of propeller and 
type of construction is determined by three principal 
factors—steady stress limits, vibration level, and pitch 
changing requirements. The first two are closely re- 
lated since any reduction in the vibration level makes 
possible an increase in the steady stress limit. Hence, 
the absence in the gas turbine of the torque impulses 
and the reciprocating forces for which piston-engine 
propellers must be designed, will permit higher steady 
stresses and make possible a reduction in propeller 


weight. However, major improvements in piston en- 


gines in the last few years, such as dynamic dampers 


and second order balancers, have made it possible to 
reduce weights per horsepower of propellers for these 
engines, and savings as great as would have been pos- 
sible several years ago are not likely. A further factor 
is that aerodynamic excitation is not reduced by the 
type of power plant. Hence, the limit on weight re- 
duction would appear to be established by the amount 
of aerodynamic excitation that must be designed for. 
At the present time a weight reduction of about 15 per 
cent appears possible. 

Effects of vibration of the propeller on the airplane 
must be considered. The development of modern 
flexible or dynamic mounts for piston engines has pro- 
gressed to the point where they are highly effective, 
but some thought has been given to eliminating this 
type of mount for turbine-propeller installations. Al- 
though propellers can be balanced statically and aero- 
dynamically to extremely close tolerances, it is prob- 
able that demands for vibrationless cockpits and cabins 
will in the future require that turbine mounts be used 
which will have sufficient flexibility to prevent trans- 
“mission of propeller order vibration. 

One other design factor that should be considered 
early in the design of turbine propellers is the type of 
deicing system. The use of heated air through a hol- 
low blade appears the best solution to the deicing prob- 
lem. The source of heated air can be from the turbine 
exhaust, the compressor discharge, or by a separate gas 
burner. Temperatures of 350-400°F. at the propeller 
appear desirable and the quantity of air required is 
estimated at 0.3 Ib. per sec. for a 15-ft. diameter three- 
blade 150 activity factor propeller having partitioned 
blades. For other sizes of propellers the airflow re- 
quirement can be scaled up or down as a first approxi- 
mation by the square of the diameter and the first 
power of the activity factor. 


TURBINE-PROPELLER CONTROL 


There are several characteristics of the turbine which 
make the control system different from that used with 
the conventional constant-speed propeller for a piston 
engine. These are as follows: 

(a) Allowable overspeed is exceedingly small, de- 
pending on the turbine design, but limits as low as | 
per cent have been proposed. 

(b) Since for maximum performance and efficiency 
it is necessary to operate the turbine wheel at the high- 
est possible temperature consistent with acceptable 
service life the maximum temperature must be held 
within close limits. 

(c) Because of the relatively high ratio of kinetic 
energy in the rotating parts to the torque available for 
acceleration and deceleration the response of the sys- 
tem to changes in load and power setting is o1ly a frac- 
tion as fast as in the piston engine. 


(Continued on page 468) 
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Tensile Properties Affecting the Formability 
of Aluminum-Alloy Sheet at Elevated 


Temperatures 


A. E. FLANIGAN,* L. F. TEDSEN,? ano J. E. DORNt 


University of California 


ABSTRACT 

The tensile properties thought to influence formability have 
been determined for 15 aluminum-alloy sheet materials tested at 
elevated temperatures. Complete load-extension diagrams have 
been obtained, and the effects of temperature, strain rate, and 
exposure time have been noted. An appreciable increase in the 
elongation over a zero gage length is found at elevated tempera- 
tures. This suggests that forming operations in which this prop- 
erty is a criterion may be facilitated at such temperatures. A 
decrease in the limit of uniform elongation is noted at elevated 
temperatures, suggesting that operations in which this property 
is a criterion may be less successful. A fiftyfold change in strain 
rate has an appreciable effect on properties at the higher tempera- 
tures, an increase in strain rate appearing to be equivalent to a 
decrease in temperature. In the case of materials subject to 
precipitation hardening, the time at temperature appears to be a 
minor factor for times ranging from 5 to 20 min. and tempera- 
tures up to at least 450°F. In certain cases unusual load-exten- 
sion diagrams were noted. 


INTRODUCTION 


In the aircraft industry aluminum-alloy sheet is often 
subjected to plastic forming, generally at room tempera- 
ture. It has been found that some of the recent high- 
strength materials exhibit poor forming characteristics 
under these conditions and it has been suggested that 
formability may be improved at elevated temperatures. 
Previous investigations! ?: * dealing with the elevated 
temperature properties of aluminum alloys have not 
included some of the recent materials and have not been 
concerned with formability. Consequently, an inves- 
tigation was undertaken to evaluate over a considerable 
range of temperatures the tensile properties related to 
formability. The program comprising a portion of the 
“restricted” Project NRC-548, was financed by the 
Office of Production Research and Development of the 
W.P.B. and was supervised by the War Metallurgy 
Committee. This paper, which presents some of the 
results, has been released for publication by the latter 
agencies. 


TENSILE DaTa, FLOW CURVES AND FORMABILITY 


The load deformation data derived in a tensile test 
are of interest in connection with formability since they 
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describe the strain hardening characteristics of the ma- 
terial and the conditions for necking and fracture. It 
is often desirable, however, to express the data in terms 
of a true stress-true strain relation (flow curve) in pref- 
erence to the more familiar but less fundamental engi- 
neering stress-strain diagram.‘ This is possible on the 
assumption that the flow curve approximates a straight 
line between the points of maximum load and fracture 
and that the incidence of maximum load corresponds 
to the beginning of necking (the limit of uniform elonga- 
tion). True stress and true strain may then be ex- 
pressed as:** 


o = s(1 + E), = log, (1 + £) 


up to and including the limit of uniform elongation, and 
subsequently as 


a = s (Ao/A), ¢ = log, (Ao/A) 


where 
o = true stress (load divided by the corresponding 
momentary cross-sectional area) ft 
s = engineering stress (load divided by the original 


cross-sectional area) 
E = engineering strain = (J — h)/l 
@ = true strain 


®Ao 


= dA/A = f dl/l up to the limit of uniform 
elongation 
= fra /A beyond the limit of uniform elonga- 
tion 


Ao = original cross-sectional area 


A = momentary cross-sectional area 
l, = original gage length 
1 = momentary gage length 


From this point of view, two of the most significant 
features of the load-extension diagram are the points of 
maximum load and fracture, and associated with these 
points are the following desirable data: 


** These relations and other aspects of the flow curve are dis- 
cussed in references 4 and 5. 
tt Here we neglect the undoubted presence of a nonuniform 
stress distribution through the cross section of a tensile specimen 
in the necked region. 
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(1) Ultimate tensile stress (defining the limit of uni- 
form elongation). 

(2) Elongation at maximum load (the limit of uni- 
form elongation). 

(3) Reduction in area at the fracture. 

(4) Load or stress at fracture. 

In the present investigation complete load-extension 
diagrams were obtained for a large number of materials 
and conditions, and the data listed above under points 
1, 2, and 3 were derived and tabulated. Considerable 
difficulty was encountered, however, in the precise de- 
termination of the load at fracture—a matter that will 
be discussed presently. 

It is not the purpose of this paper to prescribe meth- 
ods of applying tensile data in the prediction of limits of 
formability. In certain types of forming operations 
tensile data are directly applicable, whereas in other 
operations the method of application is far from simple 
and such factors as combined stresses and criteria for 
rupture and necking may be involved. 


CLASSIFICATION OF MATERIALS 


The materials investigated may be grouped as fol- 
lows according to the anticipated effects of elevated 
temperature exposure on subsequent room temperature 
properties: 

Group 1A.—Nonprecipitation hardenable alloys in 
the annealed condition. These materials may be 
formed at any temperature up to the temperatures of 
grain growth or hot shortness without significant ef- 
fect on the final properties. 

Group 1B.—Precipitation hardenable alloys in the 
annealed condition. These materials respond in a man- 
ner similar to that characterizing Group 1A. After 
forming, however they may be heat-treated to obtain 
marked improvement over the original properties. 

Group 2A.—Precipitation hardenable alloys pre- 
viously subjected to solution heat-treatment followed 
by aging at room temperature. These materials are 
subject to additional precipitation at elevated tempera- 
tures dependent upon the temperature and time of ex- 
posure. Brief exposures at temperatures below ap- 
proximately 500°F. may be expected to increase the 
room temperature yield stress. Overaging attended 
by a loss of strength may be encountered as a result of 
longer exposures and higher temperatures. 

Group 2B.—Precipitation hardenable alloys pre- 
viously subjected to solution heat-treatment followed 
by aging at elevated temperatures. Additional ex- 
posure at elevated temperatures may be expected to 
induce little additional strengthening. As in the case 
of the 2A materials, however, overaging may result at 
the higher temperatures. 


SCOPE OF THE INVESTIGATION 


The following materials were studied. All were ob- 
tained from the producers’ commercial stock. Unless 
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TABLE 1 
NOMINAL COMPOSITIONS OF THE ALLOYS TESTED * 


2S Aluminum and normal impurities. No alloying addi- 
tions 

38 1.2 per cent manganese 

52S 2.5 per cent magnesium, ().25 per cent chromium 

248 +.5 per cent copper, 0.6 per cent manganese, 1.5 per 
cent magnesium 

61S 0.25 per cent copper, 0.6 per cent silicon, 1.0 per cent 
magnesium, 0.25 per cent chromium 

R301 4.5 per cent copper, 1.0 per cent silicon, 0.8 per cent 


manganese, 0.4 per cent magnesium 

No information has been released on the precise com- 
position. The principal hardening elements are 
zinc and magnesium with smaller additions of cop- 
per and other elements 


XB75S 


* In each case the remainder consists of aluminum and normal 
impurities. 


TABLE 2 


Standard Heat-Treatments Used in Production of the Age 
Hardened Materials 


248S-T 24S alloy solution heat-treated at 910-930°F. fol- 
lowed by quenching in cold water and aging at 
room temperature 


24S-RT 24S-T strain hardened by rolling to a reduction of 
about 5.5% 

61S-W 61S alloy solution heat-treated at 960-980°F. fol 
lowed by quenching in cold water and aging at 
room temperature 

24S-T84 24S-T stretched 3!/.-4 per cent and aged 6 hours 
at 375°F. or 9 hours at 365°F. 

24S-T86 24S-RT aged 5!/2 hours at 375°F. or 8'/» hours at 
365 °F. 

61S-T 61S-W aged 8 hours at 350 °F. or 18 hours at 320°F. 


XB75S alloy solution heat-treated at 860-930°F 
followed by quenching in cold water and aging 
24 hours at 250°F. 

R301 alloy solution heat-treated at 930-950°F 
followed by quenching in cold water and aging () 
hours at 350°F. or 18 hours at 320°F 


XB75S-T 


R301-T 


otherwise indicated, all were in the bare (unclad) con- 


dition. 

Group 1A Group 1B Group 2A Group 2B 
25-0 245-0) » 24S-T 24S-T-84 
3S-0 61S—0 24S-RT 24S-T-86 

52S-0 XB75S—0 Alclad* 61S-W 61S-T 

R301-0 (clad) XB75S-T* 


R301-—T (clad 


The nominal chemical compositions of the alloys and 
the standard heat-treatments for the precipitation 


hardened materials are listed in Tables 1 and 2. In| 


each case elevated temperature tensile tests were made 
on 0.064-in. sheet to determine: 


(1) The complete load extensiom diagram 

(2) Ultimate tensile stress. 

(3) Limit of uniform elongation. 

(4). Local elongation at the fracture (based on re- 
duction in area). 

(5) Elongation in 2 in. 

(6) Effect of gage length on elongation.? 


* The alloy formerly designated X B75S by the Aluminum Com- 
pany of America is now known as 75S. 

+ The extensive data relating per cent elongation to gage length 
have been omitted from this paper. They are included in refer- 


ence 5. 
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CROSSHEAD SPEED , 


XB75S-0 
400 ‘F 


LOAD , Pounds 


4 ~< 


EXTENSION, in per 2 in EXTENSION , in per 2 in 


TIMING RATE: 


100 pulses /min 30 pulses /min 


CROSSHEAD SPEED, 53 


pounds 


xXB75S-0 
XxB75S-O 608 400 °F 


LOAD, 


70 °F 


EXTENSION, in. per 2 in EXTENSION, in. per 2 in 
TIMING RATE: 


120 pulses /sec 30 pulses /sec 
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in. / minute 


XB75S-O 658 
900 °F 
~~. 


EXTENSION, in. per 2 in EXTENSION , in. per 2 


27 pulses Amin 5S pulses /min 


in. /minute 


EXTENSION, in. per 2 in EXTENSION , in. per 2 in 


15 pulses /sec 20 pulses /sec 


Fic. 1. Typical load-extension records for XB75S-O Alclad sheet 


With the exception of items 5 and 6, each of 
these characteristics is thought to influence forma- 
bility 

Specimens were stressed transverse to the direction 
of rolling. The testing variables included temperature, 
time at temperature, and crosshead speed, the following 
values of which were employed in various combina- 
tions: 

Temperature 

Group 1: Room temperature,* 300°, 400°, 500°, 
700°, 800°, 900°F. 

Group 2: Room temperature,* 300°, 350°, 400°, 
450°, 500°F. 

Time at Temperature (including a maximum of 2 

in. for heating) 

10 min. 
5, 10, and 20 min. 


Group |: 

Group 2: 
Crosshead Speed 

3 in. per min. and 153 in. per min. 

A maximum temperature of 500°F. was used with 
Group 2 on the assumption that objectionable changes 
in room temperature properties might result from ex- 
posure to higher temperatures. The 10-min. exposure 
period was considered generally suitable for elevated 
temperature forming. In the case of the precipitation 
hardened materials of Group 2, however, additional 
spot tests were made after 5 and 20 min. The cross- 
head speeds employed bracket the range of speeds used 


in many forming operations. 


LoAD-EXTENSION RECORDS 


Load-extension diagrams for each test were obtained 
by means of a specially developed autographic method 
described in reference 6. According to this method, 
two voltages are regulated in such a manner that the 
first is proportional to the load on the specimen and 
the second is proportional to the extension in a 2-in. gage 
length. 

The resulting voltages are applied across the de- 
flecting plates of cathode-ray oscilloscope, causing 
the beam of the latter to trace out a load-extension curve 
which is recorded photographically. The beam is ex- 
tinguished at regular intervals during the course of a 
test to produce a dashed curve from which it is possible 
to determine variations in strain rate. Upon final 
parting of the specimen the beam is automatically ex- 
tinguished permanently. Typical load-extension re- 
ords are shown in Fig. 1. 

Values of maximum load and elongation at fracture 
indicated by the records were checked with values 
determined independently and were found in good 
agreement. Reproducibility was excellent. As noted 
in reference 6, however, difficulty was met in de- 
termining the precise loads at fracture. The problem 
arises from the fact that the apparatus continues to re- 
cord during the period of crack propagation following 
the initiation of failure. 


* Approximately 75°F 
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FORMABILITY OF 


STRAINING OF SPECIMEN 


Specimen dimensions are shown in Fig. 2, the design 
conforming to A.S.T.M. Standard E8-42. All speci- 
mens of each material were obtained from a single sheet. 
A specially constructed screw-driven testing machine 
affording a wide range of constant crosshead speeds 
was used at speeds of 3 in. and 153 in. per min. 
The problem of acceleration at the higher speeds was 
minimized by allowing the driving motor and the trans- 
mission system to attain speed before engaging an upper 
spherical seat in a series with the specimen. 

In spite of the constant crosshead speed, the rate of 
strain in the 2-in. gage length varied during the course 
of a test, the rate in the elastic portion of the extension 
being considerably less than that attained in the plastic 
region. In the latter portion the average strain rate 
was approximately 0.015 per sec. at the lower crosshead 
speed and 0.85 per sec. at the higher speed. The varia- 
tion during a test is attributed primarily to specimen de- 
sign. In the early stages of a test both the reduced sec- 
tion and the much longer grip section were strained 
elastically so that only a small part of the crosshead 
motion was concentrated in the gage length. Because 
of its smaller section, however, the latter flowed plas- 
tically at loads at which the remainder of the specimen 
was well below the elastic limit. Thus, an increasing 
proportion of the crosshead motion was concentrated 
in the gage length in the latter stages of a test. 


HEATING OF SPECIMENS 


Specimens were brought to the testing temperatures® 
in less than 2 min. and were maintained throughout the 
holding and testing periods with a temporal variation 
no greater than +3°F. at all temperatures employed. 
Spatial variations along the gage length were found to 
be +2°F. at 300°F. and +4°F. at 800°F. 

Two separate sources of heat were applied. The first 
was a split cylindrical furnace with three banks of heat- 
ing elements arranged vertically (Fig. 3). To secure 
rapid heating of the specimen, an auxiliary source of 
heat was provided by passing a controlled current 
through the specimen itself. This current, which was 
introduced through flexible copper leads attached to the 
grips, was gradually reduced as the furnace approached 
the desired temperature. Since the controlling pyrom- 
eter governing the furnace was set for a temperature 
slightly lower than that desired in the specimen, the use 
of a small auxiliary current allowed precise control. 
The current was reduced to zero just prior to testing, 
but due to the short duration of a test the resulting 
drop in temperature was negligible. 

Temperature measurement was made by means of 
30-gage iron-constantan thermocouples secured to the 
surface at the center of the gage length of each specimen 
by several turns of flexible glass sheathing. Such 
couples do not indicate true specimen temperature 
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Fic. 8. Elevated temperature stress-strain curves for 61S-O bare 
sheet. 


since they contact the glass wrapping, but suitable cor 
rections were established at the various temperatures. 
The corrections were small, that at 800°F. being —3°F. 


EXPERIMENTAL RESULTS 


Experimental results are shown in Figs. 4 to 33. Be- 
cause of the large number of tests involved, individual 
load-extension records have not been reproduced in this 
paper. Instead, the curves from a series of representa- 
tive photographs have been replotted for presentation 
in Figs. 4 to 18. In replotting, the timing impulses of 
the original records have been omitted. As noted 
previously, difficulty was encountered in deter- 
mining the loads at fracture from study of the rec- 
ords. As a result, there may be considerable uncer- 
tainty in some cases concerning the validity of the frac 
ture stresses shown in the replotted curves. 

Figs. 19 to 33 show the effects of testing temperature, 
crosshead speed, and, in certain cases, time at tempera- 
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Fic. 19. Effects of temperature and speed on the tensile 
properties of aluminum-alloy sheet—2S-O bare. Sheet thick- 
ness: 0.064 in.; direction: cross-grain. 


© Crosshead speed 3 in. per min. following 10 min. at tem- 
perature. Z 

@ Crosshead speed 153 in. per min. following 10 min. at tem- 
perature. 


(Each point represents average of two or more tests.) 
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Fic. 20. 
Sheet thick- 


properties of aluminum-alloy sheet—3S-O bare. 
ness: 0.064 in.; direction: cross-grain. 

© Crosshead speed 3 in. per min. following 10 min. at tem- 
perature. 

@ Crosshead speed 153 in. per min. following 10 min. at tem- 
perature. 

(Each point represents average of two or more tests. 





















































© * 
rf . 2000) 
Lresco y 
a 2% = 
Be TTS 4* 
: iam ee oi 
Ty ry 7 " 
| see-0 Sar c b 














i] 








= —+—_—_+——. ——T 
[ Lt | we 
200 400 600 800 1000 
TEMPERATURE — °F 


Effects of temperature and speed on the tensile 
Sheet thick- 




















pBs ss 88s 

XN ot 
| 

UNIFORM ELONGATION, % 

















% ELONGATION IN O2 WH. 





° 





Fic. 21. 
properties of aluminum-alloy sheet—52S-O bare 
ness: 0.064 in.; direction: cross-grain 

© Crosshead speed 3 in. per min. following 10 min. at tem- 
perature. 

@ Crosshead speed 153 in. per min. following 10 min. at tem- 
perature. 

(Each point represents average of two or more tests 
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Fic. 22. Effects of temperature and speed on the tensile 
properties of aluminum-alloy sheet—24S-O bare. Sheet thick- 
ness: 0.064 in.; direction: cross-grain. 

O Crosshead speed 3 in. per min. following 10 min. at tem- 
perature. 

@ Crosshead speed 153 in. per min. following 10 min. at tem- 
perature. 

(Each point represents average of two or more tests.) 
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Fic. 23. 
properties of aluminum-alloy sheet—61S-O bare. 
ness: 0.064 in.; direction: cross-grain. 
© Crosshead speed 3 in. per min. following 10 min. at tem- 
perature. 
@ Crosshead speed 153 in. per min. following 10 min. at tem- 
perature. 
(Each point represents average of two or more tests.) 
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Fic. 24. Effects of temperature and speed on the tensile 
properties of aluminum-alloy sheet—R301-O clad. Sheet thick- 
ness: 0.064 in.; direction: cross-grain. 

© Crosshead speed 3 in. per min. following 10 min. at tem- 
perature. 

@ Crosshead speed 3 in. per min. following 10 min. at tem- 
perature. 

(Each point represents average of two or more tests.) 
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Fic. 25. Effects of temperature and speed on the tensile 


properties of aluminum-alloy sheet—XB75S-O Alclad. Sheet 
thickness: 0.064 in.; direction: cross-grain. 

© Crosshead speed 3 in. per min. following 10 min. at tem- 
perature. 

@ Crosshead speed 3 in. per min. following 10 min. at tem- 
perature. 

(Each point represents average of two or more tests.) 
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Fic. 26. Effects of temperature and speed on the tensile 
properties of aluminum-alloy sheet—24S-RT bare. Sheet thick- 
ness: 0.064 in.; direction: cross-grain. 






































A Crosshead speed 3 in. per min. following 5 min. at tem- 
perature. 
© Crosshead speed 3 in. per min. following 10 min. at tem- 
perature. 
(_] Crosshead speed 3 in. per min. following 20 min. at tem- 
perature. 
@ Crosshead speed 153 in. per min. following 10 min. at tem- 
perature. 
(Each point represents average of two or more tests.) 
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Fic. 27. Effects of temperature and speed on the tensile 
properties of aluminum-alloy sheet—24S-T bare. Sheet thick- 
ness: 0.064 in.; direction: cross-grain. 

A Crosshead speed 3 in. per min. following 5 min. at tem- 
perature. 

© Crosshead speed 3 in. per min. following 10 min. at tem- 
perature. 

(] Crosshead speed 3 in. per min. following 20 min. at tem- 
perature. 

@ Crosshead speed 153 in. per min. following 10 min. at tem- 
perature. 

(Each point represents average of two or more tests.) 
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Fic. 28. Effects of temperature and speed on the tensile 
properties of aluminum-alloy sheet—61S-W bare. Sheet thick- 
ness: 0.064 in.; direction: cross-grain. 

A Crosshead speed 3 in. per min. following 5 min. at tem- 
perature. 

O Crosshead speed 3 in. per min. following 10 min. at tem- 

ture. 

(J Crosshead speed 3 in. per min. following 20 min. at tem- 
perature. 

@ Crosshead speed 153 in. per min. following 10 min. at tem- 


perature. 
(Each point represents average of two or more tests.) 
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Fic. 29. Effects of temperature and speed on the tensile 
properties of aluminum-alloy sheet—24S-T84 bare. Sheet thick- 
ness: 0.064 in.; direction: cross-grain. 

A Crosshead speed 3 in. per min. following 5 min. at tem- 
perature. 

O Crosshead speed 3 in. per min. following 10 min. at tem- 
perature. 

(] Crosshead speed 3 in. per min. following 20 min. at tem- 
perature. 

@ Crosshead speed 153 in. per min. following 10 min. at tem- 
perature. 

(Each point represents average of two or more tests.) 
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Fic. 30. Effects of temperature and speed on the tensile 
properties of aluminum-alloy sheet—24S-T86 bare. Sheet thick- 
ness: 0.064 in.; direction: cross-grain. 

A Crosshead speed 3 in. per min. following 5 min. at tem- 
perature. 

© Crosshead speed 3 in. per min. following 10 min. at tem- 
peraturé. 

(J Crosshead speed 3 in. per min. following 20 min. at tem- 
perature. 

@ Crosshead speed 153 in. per min. following 10 min. at tem- 
perature. 

(Each point represents average of two or more tests.) 
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Fic. 31. 
properties of aluminum-alloy sheet—61S-T bare. 
ness: 0.064 in.; direction: cross-grain. 

A Crosshead speed 3 in. per min. following 5 min. at tem- 
perature. 

© Crosshead speed 3 in. per min. following 10 min. at tem- 
perature. 

(CD Crosshead speed 3 in. per min. following 20 min. at tem- 
perature. 

@ Crosshead speed 153 in. per min. following 10 min. at tem- 
perature. 

(Each point represents average of two or more tests.) 
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Fic. 32. Effects of temperature and speed on the tensile 
properties of aluminum-alloy sheet—R301-T clad. Sheet thick- 
ness: 0.064 in.; direction: cross-grain. 

A Crosshead speed 3 in. per min. following 5 min. at tem- 
perature. 

© Crosshead speed 3 in. per min. following 10 min. at tem- 
perature. 

(J Crosshead speed 3 in. per min. following 20 min. at tem- 
perature. 

@ Crosshead speed 153 in. per min. following 10 min. at tem- 
perature. 

(Each point represents average of two or more tests.) 
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Fic. 33. Effects of temperature and speed on the tensile 
properties of aluminum-alloy sheet—XB75S-T bare. Sheet 
thickness: 0.064 in.; direction: cross-grain. 
A Crosshead speed 3 in. per min. following 5 min. at tem- 
perature. 


© Crosshead speed 3 in. per min. following 10 min. at tem- 
perature. 

(_] Crosshead speed 3 in. per min. following 20 min. at tem- 
perature. 

@ Crosshead speed 153 in. per min. following 10 min. at tem- 
perature. 


(Each point represents average of two or more tests.) 
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Fic. 34. Typical load-extension record for 52S-O sheet at room 


temperature. 


ture on the various tensile properties. In general, each 

point represents an average value based on two or more 

tests. The data were obtained in the following man- 

ner: 

Ultimate Tensile Stress-—From the maximum load 
shown by the'load-extension record. 

Limit of Uniform Elongation.—From the strain corre- 
sponding to maximum'load on the load-extension record. 
Both plastic and elastic components of strain are in- 
cluded. The results are plotted in terms of true strain 
as well as engineering strain. 

Elongation in 2 \In.-From measurement of the frac- 
tured specimen. ~~ “ph 

Local Ductility —(equivalent to the elongation over 
a zero gage length at the fracture)—-From measurement 
of the reduction in area at the fracture. Measurements 
were made on both halves of each specimen by means of 

‘a toolmaker’s microscope and were found to agree 
closely. The results are plotted in terms of local duc- 
tility and true plastic strain according to the relations. 


‘ 


Local ductility in per cent = 100[(Ao0/A) — 1] 
True plastic strain at fracture = log, (Ao/A) 


Because of the uncertainties concerning the load at 
fracture, fracture stresses have not been included in 
Figs. 19 to 33. 


Effect of Temperature on Tensile Properties 

Figs. 19 to 33 indicate that in all materials relatively 
little change in tensile properties occurs at temperatures 
up to 300°F. At higher temperatures the changes are 
more pronounced, and in some cases hot shortness is 
evident at 800°F. and 900°F. In general, the effects 
of temperature are those anticipated for the various 
groups of materials and, therefore, no further discussion 
of this point will be made. 

It is of interest to note that the elongation in 2 in. 
may be a poor criterion of formability. This is not un- 
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expected since, in reality, it is not a fundamental 
property but is an arbitrary determination falling be- 
tween values of the more fundamental elongation over 
a zero gage length and the limit of uniform elongation, 
These properties vary considerably with an increase 
in temperature, the former tending to increase and the 
latter to decrease. It may seem from Figs. 32 and 33 
that in the cases of R301-T and XB75S-T there is a 
considerable increase in the elongation over a zero gage 
length at temperatures above 400°F., whereas the in- 
crease in the elongation in 2 in. is much smaller. In the 
case of the Group 2A materials (Figs. 26 to 28) it is 
seen that at the higher temperatures there may be an 
appreciable increase in local elongation even when the 
elongation in 2 in. exhibits a decrease. This implies 
that in cases where local ductility:is a criterion the form- 


-ability at elevated temperatures may improve much 


more than indicated by values of the 2-in. elongation. 
Conversely, in cases where the limit of uniform elonga- 
tion serves as criterion, elevated temperature formabil- 
ity may be inferior to that indicated by the elongation 


in 2 in. 
Effect of Time at Temperature 


In the cases of 24S-T and R301-T (Groups 2A and 
2B, respectively), tests involving 5- and 20-min. holding 
periods were made in addition to tests with the standard 
time of 10 min. It will be noted in Figs. 27 and 32 
that at all temperatures considered the effect of time is 
rather small. The remaining materials of Group 2 were 
tested after 5- and 20-min. exposures at a single tempera- 
ture, usually 450°F. In these cases, again, the effect 
of time is small. It is concluded, therefore, with regard 
to exposure time, that the elevated temperature proper- 
ties listed in this paper may be used with reasonable as- 
surance in commercial forming operations. 


Effect of Strain Rate 


It is of interest to note the effect of strain rate on the 
shape of the stress-strain curve. In general, it may be 
seen that the curves derived at the higher speed at a 
given temperature tend to resemble those derived at 
the lower speed at a lower temperature. This tendency 
is especially noticeable in the case of Group 1 materials 
where precipitation effects are not encountered. 

Recently Zener and Holloman have proposed an 
equivalence of the effects of changes in strain rate and 
in temperature upon the stress-strain relation.‘ Al- 
though the range of rates employed in the present 
tests is small, it is possible that an examination of the 
data may be used to check the proposed equivalence. 

The effect of strain rate on the individual tensile 
properties is shown in Figs. 19 to 33. Again it is seen 
that the high-speed properties at a given temperature 
correspond in general to low-speed results obtained at a 
lower temperature. In addition, it is apparent that 
the effect of speed is rather small in all materials at 
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FORMABILITY OF 


room temperature and at 300°F. but becomes increas- 
ingly important at higher temperatures. 


Unusuat Features of the Load-Extension Records 


Several interesting features of the original load-ex- 
tension records are noted in the following: 

(1) A markedly pronounced knee is observed in the 
diagram of 52S-O at room temperature (Fig. 34). In 
addition, it may be seen that the entire curve consists 
of a series of such knees. During testing it was ob- 
served in this case that elongation was accomplished 
through a series of local yieldings. The same type of 
load-extension diagram is found to a lesser extent with 
24S-O, R301-O, and XB75S-O (see Fig. 1). 

(2) Inthe case of the four materials just mentioned, 
sharp peaks are found on the load-extension diagrams 
derived at the higher temperatures. Apparently there 
is a connection between this effect and the appearance 
of knees in the diagrams of the same materials at room 
temperature. The peak is illustrated for 52S5-O in 
Fig. 35 (See also Fig. 1.) 

(3) Fig. 36 typifies the appearance of diagrams de- 
rived from the testing of 24S-RT at room temperature 
at the lower crosshead speed. In each of the four tests 
a ‘“‘screeching’’ noise was heard just before fracture. 
The curves exhibit an unusual appearance in this re- 
gion. Neither the ‘“‘screeching’’ nor this type of dia- 
gram was noted in the other materials. 


’ 


SUMMARY AND CONCLUSIONS 


(1) Tensile properties thought to influence form- 
ability have been determined for 15 aluminum-alloy 
sheet materials tested at elevated temperatures. Com- 
plete load-extension diagrams have been obtained for 
each condition and the effects of temperature, strain 
rate, and exposure time have been noted. 

(2) An appreciable increase in the elongation over a 
zero gage length is found in all materials at elevated 


temperatures. This suggests that forming operations 





Fic. 35. Typical load-extension record for 52S-O sheet at 800 °F. 
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Fic. 36. Typical load-extension record for 24S-RT sheet at room 


temperature 


in which this property is a criterion may be facilitated 
at such temperatures. 

(3) A general decrease in the limit of uniform elon- 
gation is noted in all materials at elevated tempera- 
tures. 
in which this property is a criterion may be less success- 
ful. 

(4) A fiftyfold change in strain rate is found to have 
an appreciable effect on properties at the higher tem- 


It may be expected, therefore, that operation 


peratures. The rate of deformation may, therefore, 
be expected to be a factor in formability. 

(5) In general, the effect of an increase in strain 
rate is equivalent to that of a decrease in temperature. 

(6) In the case of materials subject to precipitation 
hardening, the time at temperature is found to be a 
minor factor for times ranging from 5 to 20 min. This 
is true in all cases for temperatures up to at least 
450°F. 

(3a 
grams were noted. 


several cases unusual load-exteusion dia- 


These have been described. 
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Propellers for Aircraft Gas Turbines 


(Continued from page 456) 


(d) To start the turbine, the r.p.m. has to be 
brought up to a fairly high speed (20-25 per cent). 
To keep the power required for starting to as low a value 
as possible, some means of putting the propeller in flat 
pitch is desired. 


(e) In the landing approach it is necessary that 
power be available almost instantaneously, at least in 1 
sec. This requires that the turbine operate at near its 
rated speed. Control of the power plant to provide 
zero thrust or some preferred negative thrust appears 
necessary. 

No attempt is made in this paper to go into the vari- 
ous possible types of control systems to meet the above 
requirements. It is merely stated that it appears that 
a different type of control system than the constant- 
speed governor system used on piston engines will be 
necessary. One of the most promising arrangements is 
to have the governor on the fuel regulator and to con- 
trol the propeller fromsome indirect measure of tempera- 
ture, such as fuel-air ratio or torque to change pitch to 
absorb the desired amount of power. When a direct 


method of turbine temperature measurement is evolved 
which is sufficiently rapid, accurate, and reliable for 
control purposes, this can be used directly to provide 
the basic propeller control. 


CONCLUSION 


This paper has discussed the general problems as- 
sociated with the design of propellers for aircraft gas 
turbines and the application of turbine-propeller 
power plants to high-speed airplanes. The conclusion 
is reached that at speeds of 600 to 700 m.p.h. propeller 
efficiencies of over 70 per cent appear possible, particu- 
larly with unconventional configurations, such as 
swept-back blades and small-diameter multiblade pro- 
pellers. This high-speed performance is combined with 
superior performance for take-off, climb, and cruising. 
The characteristics of propellers for use on long-range 
aircraft having speeds of 400-500 m.p.h. are also dis- 
cussed and it is shown that to obtain best cruising 
efficiency and good take-off characteristics, propellers 
of high solidity appear desirable. 
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Flight-Test Investigation of the Strength of 
_a Square-Tipped Wing 


LEONID A. MINOR* 


Grumman Aircraft Engineering Corporation 


SUMMARY 


Loads measured in the outer portion of a square-tipped wing 
during flight differ substantially from the design loads calculated 
according to ANC-1 specifications. A rapid increase of loads per 
gat the tip is observed at certain speeds and accelerations because 
of stall or compressibility effects. In bending, the critical wing 
section shifts from the tip, at high altitude where the stall and 
buffet boundaries can be attained to the root, at low altitudes 
where, in general, the operating restrictions limit the normal 
accelerations well below these boundaries. 

Test data will be presented and analyzed, and conclusions will 
be drawn relative to the strength that is required at different wing 
stations in order to avoid critical failures if design limit accelera- 
tions are exceeded in flight. 


INTRODUCTION 


Ws STRUCTURES are generally designed on the 
basis of calculated span load distributions for a 
number of specified flight conditions. Symmetrical 
flight conditions usually set minimum strength require- 
ments for a constant design load factor throughout a 
specified speed range. Lack of agreement between cal- 
culated design loads and actual flight loads may not 
necessarily cause a failure in flight below the design 
load factor if sufficiently high design margins of safety 
are provided throughout the wing span to cover the 
probable inaccuracies of the applied load calculations. 
However, arbitrary_margins of safety result automati- 
cally in a waste of structural weight. 

On a combat aircraft one cannot expect rigid g re- 
strictions to be scrupulously observed and high positive 
design margins do not necessarily increase the safety of 
the airplane in proportion to the weight expenditure in- 
volved, unless they are distributed along the wing span 
in such a manner that, should a failure occur, it will not 
cause the total loss of the airplane. The same considera- 
tions hold true for commercial aircraft on which de- 
sign load factors can be accidentally exceeded. This 
may occur, for example, under severe gust condi- 
tions. 


The desirability of achieving a localized failure led the - 


Grumman Aircraft Engineering Corporation to experi- 
ment with wing-tip safety joints for the purpose of 
demonstrating their usefulness as a safety device. The 
wing of a production airplane was modified to break in 
symmetrical flight just outboard of the middle aileron 

Presented at the Design Session, Fourteenth Annual Meeting, 
1.A.S., New York, January 29-31, 1946. 

* Structures Flight Research Engineer. 


hinge under normal accelerations corresponding to 50 
per cent of the limit design load factor. Flight tests 
showed the practical possibilities of this design feature. 
The two wing tips and the aileron tips broke off simul- 
taneously during a symmetrical pull-out. The airplane 
was found to be perfectly maneuverable with its wings 
clipped, and a safe landing was made under adverse 
conditions (20 m.p.h. cross wind). Flights were also 
made with only one wing tip off and showed that the 
airplane in that configuration had satisfactory flying 
characteristics and could make safe take-offs and land- 
ings. 

Extensive instrumented flight tests were conducted 
thereafter by Grumman’s Structures Flight Research 
Group on a recent fighter model airplane featuring 
square-tipped wings with the purpose of establishing 
speed and normal acceleration restrictions in symmetri- 
cal flight at any altitude based on the actual strength 
of the wing. Loads were measured at three stations of 
the wing—at the root, at mid-span, and near the tip— 
in order to establish the desirable relationship of 
strength between these three sections so that the tip 
section would be critical under any symmetrical flight 
condition. The results of this investigation provide the 
basic data for the present paper. 


TEst AIRPLANE AND TEST EQUIPMENT 


The test airplane was a single-engined fighter of con- 
ventional design with square-tipped wings of aspect 
ratio 5.15 and taper ratio 2.24, with a modified NACA 
23018 section at the root and an NACA 23009 section 
at the tip. The gross weight for all pull-outs was kept 
constant within 1 per cent, while the c.g. travel was 
restricted within '/. per cent of -M.A.C. 

The airplane was instrumented with a set of strain- 
gage equipment and a twelve-channel oscillograph. 
The records were taken on an 8-in. wide paper tape, 
with some traces moving as much as 7 in. during high g 
pull-outs. Normal acceleration, elevator position, 
aileron position, stick force, and wing bending moments, 
shears, and torsions at sections situated at approxi- 
mately 12.5, 50, and 80 per cent of the semispan were 
recorded. The strain gage equipment was supple- 
mented by a motion-picture camera installed in the 
cockpit and viewing the air-speed indicator, the altim- 
eter, and the visual accelerometer on the instrument 


panel. Both sets of equipment were operated from a 
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single switch, and a timing device was used for syn- 


chronization. 


Instrument Calibration 

A highly sensitive bearn-type accelerometer cali- 
brated on a specially designed pendulum was used. The 
accuracy of the recordings was estimated to be +0.10g. 

Bending, shear, and torsion gages at each wing station 
were calibrated simultaneously. Because the gages 
were not independent, three calibration loads plus a 
check load were used for each station, and loads were 
expressed linearly in terms of trace displacements of all 
three gages. The air-speed indicator and the altimeter 
were calibrated for instrument and position error. Air 
speeds were also corrected for impact pressure com- 
pressibility error. 

All calibrations were repeated at regular intervals in 
order to check variations of sensitivity of the various 


instruments. 


FLIGHT-TEST PROGRAM 


A first series of flights was made at high altitude to 
determine the low-speed maximum lift coefficients and 
the effects of compressibility on the maneuverability of 
the airplane. 

Low-speed accelerated stalls and compressibility 
shock stall cause sudden changes in load distributions 
and severe vibrations (buffeting) which may result in 
excessive structural stresses if these effects are encoun- 
tered at high pull-out load factors. Because of the re- 
duced speed of sound at high altitudes, compressibility 
effects appear at low values of indicated air speed and 
load factor so that these effects can be investigated 
under safe conditions before proceeding to high ac- 
celeration pull-outs at low altitude. 

Subsequent tests showed that the effect of altitude 
on the speed of sound could be eliminated by expressing 
pull-out data in terms of airplane normal force coeffi- 
cient Cv, and Mach number rather than load factor and 
speed. Consequently, the load factors and speeds at 
which changes in load distributions were liable to appear 
at low altitudes could be estimated. 

Thereafter, wing loads were measured during pull- 
outs and accelerated stalls at altitudes ranging from 
5,000 to 15,000 ft. 

Records were taken for a total of nearly 200 acceler- 
ated stalls, pull-outs, and high-speed dives. 


INTERPRETATION OF FLIGHT-TEST DATA 


Aerodynamic Limitations of the Airplane 


The flight range of an airplane in symmetrical flight, 
expressed at a given altitude in terms of load factor and 
speed, is limited by the maximum attainable airplane 
normal force coefficient Cv, at any given Mach number. 
The Cw, coefficient is proportional to the normal com- 
ponent of the air load acting on the whole airplane and 
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can be readily determined from the normal acceleration 
measured at the center of gravity of the airplane and the 
air speed. 

The maximum attainable Cy, does not necessarily 
correspond to the wing Cz» :., since the ability to attain 
high lift coefficients in flight is also dependent on the 
pilot effort. The magnitude of the stick forces often 
has more significance than the theoretical Cy,,,, of the 
wing, especially when high stick forces are related to in- 
creased longitudinal stability or a loss of elevator effec- 
tiveness. 

At Mach numbers below WM = 0.55, the maximum 
attainable Cy, is governed entirely by the wing Cpp,.- . 
Sufficient longitudinal control is available to attain the 
wing CL». :, and stall the airplane. Therefore, the stall 
boundary can be readily determined at a given altitude 
in terms of attainable Cy, versus Mach number. 

Beyond a Mach number of MW = 0.55, compressibility 
effects become increasingly important. As the angle of 
attack is increasing during a pull-out, a rearward shift 
of the wing center of pressure occurs as a consequence of 
a partial compressibility stall. This shift of center of 
pressure appeared in wing-torsion measurements. A 
further increase of Cz is limited by magnified stick 
forces because of increased static longitudinal stability 
and reduced elevator effectiveness. Vibrations appear 
in the tail (buffeting) and spread gradually through the 
whole structure of the airplane. The normal accelera- 
tion may still increase slightly as the speed is rapidly 
reduced by the large compressibility increment in drag. 

The histories of such pull-outs were plotted on a Cy, 
vs. M diagram, and the envelope of all the curves thus 
obtained was defined as the buffet boundary of the air- 
plane. 

Flights up to the stall and buffet boundaries were 
made at altitudes ranging from 5,000 to 30,000 ft. and 
at Mach numbers up to / = 0.79. Agreement between 
data obtained at different altitudes was excellent, all 
points falling within the range of experimental errors 
on acceleration and air-speed measurements, thereby 
justifying the Cy, vs. ./ mode of presentation that is 
independent of altitude and gross weight. 


Presentation of Wing-Load Data 


Load coefficients—defined as wing bending moments, 
shears, and torsions divided by g = '/29 *—-were plotted 
for each maneuver against the airplane normal force 
coefficient Cy,, rather than load factor, in order to elimi- 
nate the effect of altitude on compressibility phenomena 
and to compensate for small variations of gross weight. 
Those load coefficients derived from flights made be- 
tween 5,000 and 15,000 ft. were found to be independent 
of altitude within the range of experimental errors and 
to vary only with Mach number and Cy,. 

By using data from several pull-outs made at the 
same speed and altitude, load coefficients at constant 
Mach number were plotted against Cyv,. The wing 
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ic. 1. Wing bending moment coefficients at 80 per cent of 
semispan. 


vibrated in the first symmetrical bending mode during 
stalls or after the start of buffeting. Vibratory loads 
were taken into account and added to the measured 
steady loads. 

Typical bending moment load coefficient curves are 
presented on Fig. 1. Wing loads at any given altitude 
and load factors were determined from such curves and 
plotted against equivalent air speed (Figs. 2 and 3). 


Discussion of Flight Test-Load Data 


Wing-load coefficients were plotted against the air- 
plane normal force coefficient Cy, rather than against 
the wing-lift coefficient Cz, because the former coeffi- 
cient is easily determined from m and V, without requir- 
ing delicate measurements of tail loads and because it 
may be readily reconverted into normal accelerations 
and load factors. Therefore some inconsistencies in the 
wing-load data cannot be avoided because of balancing 
and dynamic tail loads that vary according to the man- 
ner in which the maneuver is performed. The results 
were sufficiently consistent for the purpose of this in- 
vestigation as long as the pilot did not vary his tech- 
nique for pull-outs performed in a given speed range. 

It is felt that the wing-load data interpreted in this 
manner is illustrative of normal flight conditions, and 
further refinements are justified only where more accu- 
rate comparison with theory is required. 


Another particular feature affecting the wing loads 
was the symmetrical upfloat of the ailerons which in- 
variably occurred at high angles of attack. The conse- 
quent shift of the center of pressure inboard was, of 
course, not considered in the calculation of applied 
loads. For normal force coefficients lower than Cy, = 
0.80 the aileron upfloat was negligible and could be dis- 
regarded for all practical purposes. 

The increased effects: of elastic twist on the wing 
spanwise distribution were taken into account whenever 
load data obtained at a certain altitude were used to 
predict loads at lower altitudes and higher speed. 
This correction was found to be relatively small as com- 
pared to the deviation of flight-test loads from calcu- 
lated loads, since the wing of the test airplane was de- 
signed for exceptionally high rolling performance re- 
quiring great torsional stiffness. 


METHOD OF DETERMINATION OF THE WING STRENGTH 


For a wing of conventional design, the bending 
strength at any station is usually only slightly affected 
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Fic. 2. Wing bending moments at 80 per cent of semispan at 
10,000 ft. 
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by torsional loads for symmetrical flight conditions. 
A typical strength diagram representing allowable 
bending moments plotted against allowable torsional 
moments for various elements of a wing structure at a 
particular station is shown on Fig. 4. As long as the 
torsional moments remain within the range AB, the 
bending strength remains unaffected, and the strength 
of the wing in terms of flight load factors can be readily 
determined from Figs. 2 and 3 and similar plots for 
various wing sections at any given altitude and speed. 

If the bending strength of the wing is affected by tor- 
sional loads, it can still be determined by a simple pro- 
cedure. For each altitude, a grid made of two families 
of load curves—one for constant normal accelerations 
and the other for constant speeds—is drawn on a 
strength diagram similar to the one presented on Fig. 4. 
Intersection of load curves with strength lines will de- 
termine the permissible load factor at any speed and 
also indicate which item in the wing structure is critical. 
Strength lines based on analysis or static test strength 
data and flight-test loads were determined by that 
method for the three wing stations and were plotted on 
V-n diagrams for various altitudes (Figs. 5 and 6). 

At sea level, strength lines for all three stations fall 
well below the peak of the estimated buffet boundary. 
As the altitude increases, the strength lines approach 
the peak of that curve, and above a certain altitude, 
which varies for each station, no strength lines can be 
drawn. This occurs when aerodynamic limitations pre- 
vent the airplane from attaining sufficient load factors 
to produce limit loads on the structure. 


DEVIATIONS FROM CALCULATED DESIGN LOADS 


Sample calculations of applied loads for the design of 
the wing were made by standard methods neglecting 
stall and compressibility effects for symmetrical flight 
conditions. Two typical design limit load factors were 
picked with the purpose of illustrating deviations of 
flight loads from calculated loads for a high g airplane 
and for a low g airplane having the same wing loading. 
In both cases the positive low angle-of-attack design 
condition was found to be critical for the bending 
structure throughout the entire span of the wing. In- 
crements An of allowable normal accelerations, based 
on flight-test loads, over the design load factor m were 
plotted against air speed at various altitudes (Figs. 7 
and 8), with the assumption that the wing was designed 
for zero margin of safety along the entire span. 

These curves follow the same pattern for all the three 
wing stations, with only one exception, which will be 
discussed later. The allowable normal accelerations 
increase with speed from the stall boundary to a peak 
that is attained between Mach numbers of 0.50 and 
0.56 and then decrease as the buffet boundary is ap- 
proached and the wing starts vibrating. 

For low load factors the root section is more critical 
than the tip section at low altitudes. As the altitude 
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increases and as the strength lines come closer to the 
top of the buffet boundary, the tip section becomes rela- 
tively more critical than the root. This is demonstrated 
on Figs. 7 and 8, if one considers that a high-strength 
wing at a low altitude is comparable qualitatively to a 
low design strength wing at a high altitude. 

A particular feature of the test airplane, mentioned 
earlier in this paper, is the symmetrical upfloat of the 
ailerons which occurs at high lift coefficients affecting 
the spanwise load distributions. For lift coefficients 
lower than Cr = 0.8, the upfloat measured in flight 
was too small to affect the wing loads to any appreci- 
able extent. This was verified by calculations based on 
wind-tunnel data. At Mach numbers greater than M 
= (0.60, no upfloat was noticeable regardless of altitude 
and normal acceleration. Although the aileron upfloat 
affected the wing loads at high lift coefficients, general 
conclusions regarding the wing loads can still be drawn 
in the vicinity of the stall boundary. The shift of the 
center of pressure inboard due to the aileron upfloat had 
a greater effect on loads at the tip section and at mid- 
span than at the root, since the aileron area aft of the 
hinge line accounts for 24 per cent of the wing area out- 
board of the tip section against 15 per cent at mid-span 
and 7 per cent of the root. Despite the aileron upfloat, 
the bending moments per unit load factor increase at 
the tip, as well as at the root, as the stall boundary is 
approached. This indicates a considerable shift of the 
center of pressure outboard for the portion of the wing 
outboard of the tip section, which overbalances the 
shift inboard due to the aileron upfloat. However, this 
shift outboard is not sufficient to compensate the effect 
of the aileron upfloat at the mid-span station at high 
lift coefficients. ; 

Therefore, it can be presumed that the elimination 
of the aileron upfloat on the subject wing will re- 
sult in a proportionately greater increase of bend- 
ing moments at the tip than at mid-span and at the 
root in that order. The changes in the spanwise load 
distributions in the stall region will then be similar to 
those occurring near the buffet boundary. 


DESIGN CRITERIA FOR A SQUARE-TIPPED WING 


- Design Specifications and Safety Criteria 


The stress analysis of the bending structure of a wing 
is a relatively simple problem in comparison with the 
analysis of the shear structure and is generally far more 
reliable. Design margins of safety can be reduced to a 
minimum provided the design loads are known with suf- 
ficient accuracy. It is customary, therefore, to design 
wing bending structures as closely as possible, while 
relatively high margins are provided for shear struc- 
tures. The maximum saving of weight is achieved, of 
course, in the inboard portion of the wing. Near the 
tip, the structure is often overstrength because of 
production considerations, requiring the use of mini- 


10000 FEET 





An/n 


ROOT SECTION 
--- MID SPAN SECTION 
---- TIP SECTION 
2800 ee 


An/n 





EQUIVALENT AIRSPEED Ve MPH 


Fic. 7. Increments of allowable load factors over design load 
factors for a 7.0g design. 


mum gage materials that may overrule weight con- 
siderations. 

Government agencies, as a general rule, specify mini- 
mum strength requirements for the entire structure 
(which must be fulfilled for a number of arbitrary de- 
sign conditions). The location of the critical section is 
left to the contractor’s discretion, provided the failure 
does not occur below a guaranteed load factor. 

Specified strength requirements are the result of a 
compromise between such factors as maneuverability, 
operating conditions, efficiency, and durability. There- 
fore, design load factors can be exceeded in flight acci- 
dentally, and it is in the interest of both purchaser and 
the manufacturer to avoid the total loss of the airplane 
under such circumstances. The square-tipped wing 
offers distinct possibilities in this respect. 

A square-tipped wing is characterized by the nature 
of the changes in the spanwise load distribution in the 
vicinity of the stall and buffet boundaries.. The shift 
of the center of pressure outboard which occurs under 
these conditions results in a relatively greater increase 
of bending moments at the tip than at the root. For a 
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high-performance airplane, the consequence is a large 
variation of permissible normal accelerations with speed 
at a given altitude, which is quite uneconomical with 
the present type of design strength specifications. 
Straightening of strength lines, such as shown on Fig. 4, 
can hardly be achieved simultaneously at both ends of 
the speed range, since different factors contribute in each 
case to the outboard shift of the center of pressure. By 
increasing, for instance, the thickness ratio taper, some 
improvements can be obtained at low speed, although 
at the risk of a tip stall. However, this will increase the 
outboard shift of the center of pressure at high speeds, 
because the critical Mach number for thinner tip sec- 
tions will be much greater than for the thicker root sec- 
tions. 

Conversely, a reduction of the thickness ratio taper 
will cut down the shift of the center of pressure out- 
board at high speeds and increase it at low speeds. In- 
creased chord taper ratio should reduce the variations 
of allowable load factors throughout the speed range, 
but this solution is not always practical. 

In some cases, major considerations, such as low- 
speed flying characteristics of the airplane or overall 
dimension requirements, will dictate the use of square- 
tipped wings. 

The typical cases of two airplanes having the same 
wing loading and different design load factors are con- 


sidered. 


High g Design 


In the case of a high g design wing, the maximum nor- 
mal accelerations that can be attained in flight are lim- 
ited by the stall and the buffet boundaries of the air- 
plane except at low altitude in the intermediate speed 
range, where, provided that no other item of the struc- 
ture is critical, the strength of the wing is the restricting 
factor (Figs. 5 and 6). If the tip section can be made 
critical around a Mach number of 0.52 at sea level (in 
the case of the NACA 230 series airfoils), it will remain 
critical at any speed and at any altitude as evidenced 
by Figs. 4 and 5. Furthermore, when the normal ac- 
celeration increases beyond the value corresponding to 
the limit strength of the tip section, the bending mo- 
ments at that section will increase proportionately 
more than the bending moments at the mid-span or at 
the root sections, thus providing an additional safe- 
guard against a failure inboard of the wing-tip section. 
The maximum desirable strength of a picked section 
outboard of the middle aileron hinge can be determined 
from flight tests at a safe altitude in order to utilize 
fully the strength of the inboard portion of the wing. 

The design procedure would be to determine first the 
minimum strength of the tip section necessary to meet 
the specification requirements at each end of the speci- 
fied speed range for the limit design load factor, to 
determine next the permissible limit load factor for a 
tip of that strength at Mach number of about 0.52 at 
sea level, and to design the wing for this condition. 
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For such design condition the margins of safety or 
limit loads for the whole wing can be reduced to a mini- 
mum, and a small negative margin of safety tip joint 
(possibly as small as —5 per cent) will be sufficient to 
meet the safety requirements for symmetrical flight 
conditions. 


Low g Design 


In the case of a low g design wing the problem is es- 
sentially the same. However, since the strength of the 
wing in terms of g’s is well below the aerodynamic lim- 
its of the airplane at low altitude and at intermediate 
speeds, the relative strength of the tip will not decrease 
as the loads increase beyond limit as in the case of a 
high g design. Greater care is necessary in the deter- 
mination of applied loads, and comparatively weaker 
tip sections may be required. 

The design of a wing for a low g, high-speed airplane 
presents one major problem that is not stressed suffi- 
ciently on Fig. 8. If the torsional strength of the wing 
is more in line with the bending strength requirements 
than in the case of the hypothetical 4.0g wing con- 
sidered in the present analysis, an important increase of 
wing loads due to twist must be expected at high speeds. 
This will result in a considerable variation of strength 
in terms of g’s for the inboard portion of the wing 
throughout the speed range. It will also result in a 
structurally uneconomical design. Any change in wing 
plan form which will increase the area of the inboard 
portion of the wing to the detriment of the tip will pro- 
duce favorable results. Square-tipped wings should 
therefore be unadvantageous in that respect. 


4 CONCLUSION 

Square-tipped wings are subject to important span- 
wise shifts of the center of pressure with speed and lift 
coefficient and should prove structurally uneconomical 
on airplanes designed for a constant load factor through- 
out a specified speed range. However, when featured 
on a high g airplane designed to operate close to its 
aerodynamic limits, the variations in spanwise distribu- 
tion can be utilized to obtain a localized safe failure at 
the wing tip in case of accidental overload. This can 
be achieved without serious weight penalty and without 
restricting, to any appreciable extent, the operating 
range of the airplane in terms of speed and load factors. 
If sufficient aileron area is left inboard of the critical 
section, which should be located outboard of the middle 
aileron hinge, the possibility of emergency landings with 
clipped wings will greatly increase the safety of the air- 
plane. 

Calculated wing loads are unreliable in the immediate 
vicinity of the aerodynamic limits of the airplane, and 
pertinent data relative to the maximum strength of a 
predetermined tip section where the safety joint is 
located can be best derived from flight tests. 
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Supersonic Wave Drag of Thin Airfoils 


ALLEN E. PUCKETT* 
California Institute of Technology 


INTRODUCTION 


FINITE BODY moving in a nonviscous compressible 

fluid, at a speed greater than the speed of sound, 
will experience a drag, generally called “‘wave’’ drag. 
This drag may be computed in special cases by assum- 
ing that the disturbance velocities produced by the 
body are small compared to the velocity of the body. 
With this assumption of small perturbations, the equa- 
tions of motion for the fluid may be linearized, and 
solved fairly easily. This was done by Ackeret! for 
two-dimensional thin airfoils (i.e., infinite span). For 
slender symmetric bodies, the solution was given by 
von Karman and Moore.” At the suggestion of von 
Karman, the methods have been extended to the case 
of a three-dimensional body that is extremely thin in 
one dimension, such as an airfoil. 


GENERAL EQUATIONS 


It will be assumed that the motion is isentropic and 
nonviscous, with no circulation, so that the airfoil is at 
zero lift. It is also assumed that the motion is steady 
state, with free stream velocity U in the x-direction. 
The airfoil is supposed to be extremely thin in the z- 
direction. Let the local velocities in the x, y, 2 direc- 
tion, respectively, be u’, v’, w’; with the assumption 
of small perturbations, we then have 


u’= U+u4u u = 0¢/O0x 
v =v v = 0¢/dy 
w’ =w w = 0¢/dz 


where u, v, w are small compared to U. The equation 
for the potential, ¢, of the perturbation velocities then 
becomes 


oO” a) 
— M?) tote ae — ( 


where \/ = U/a, with a equal to the speed of sound in 
the free stream.’ 

The boundary condition for a thin airfoil is simply 
that the local velocity at the airfoil surface be parallel 
to the surface at any point, i.e., 


d, = w/U (on the airfoil) (2) 


where \, is the slope of the airfoil surface measured in 
the x-direction. If the airfoil is considered to be essen- 
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tially in the x-y plane, then, within the limits of this 
approximation, it is sufficient to satisfy Eq. (2) in the 
plane z = 0. 

An elementary solution to Eq. (1) is given by 


c 

oa Side eee gees eee - mer Sg (3) 
Vv (x — £&)? — Bl(y — n)? + 27] 

where 8? = M*? — 1. That this is a solution to Eq. 


(1) can be verified directly by substitution. However, 
it can also be shown to represent the potential due to a 
single source of strength proportional to c, located at 
point (é, n) in the x-y plane.‘ If the Mach Number, 
M, goes to zero, it is observed that Eq. (3) takes the 
familiar form of the source potential for an incompres- 
sible fluid. It will also be noticed that Eq. (3) becomes 
imaginary for points lying outside the cone defined by 


(y — n)? + 2? = (1/8%)(x — £)? (4) 


But this is exactly the cone with vertex at (&, »), and 
semivertex angle equal to the Mach angle, 6, where 


6 = sin— (1/M) 
(1/+/M? — 1) 


II 


(1/8). (See Fig. 1.) 


since tan 6 = 





This cone, called the Mach cone, is the instantaneous 
wave front that would be produced by a small disturb- 
ance at (£, 7), moving at Mach Number .J/ into still air. 
It is also, therefore, the stationary wave front produced 
in a flow at Mach Number J/ by a disturbance fixed at 
(&, ). 
effect of the disturbance at (£, 7) should not be felt out- 
side of the cone defined by Eq. (4). Since Eq. (4) ac- 
tually defines two cones, the disturbance will be felt 
“after- 


It is, therefore, physically reasonable that the 


only in the cone where x > £, which is called the 
cone.’’ The potential @ of Eq. (3) will be taken as 
zero everywhere outside this cone. 

A more general solution to Eq. (1) can be constructed 
by superposition of potentials of the type of Eq. (3), 
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with c a function ot £, 7. A particular form will be con- 
sidered, in which c is a surface source element, and the 
elementary potential is 
q dé dy 
V(x — 8? — By — 0)? + 2] 
Here q is the source strength per unit area. The verti- 
cal velocity, 0¢/0z, due to this source element is zero 


everywhere in the x-y plane except on the atea d§ dy; 
for, by differentiating Eq. (5), we obtain 


aw gB?z dé dn 

= (do) agli 2 - 2 anaes 
oz {(@ — 8)? — Ally — »)* + 2°} 
and at z = 0, this expression is zero. But this is valid 
only for x — & large compared to dé; in other words, 
for points exterior to the element dé dyn. For a point 
x, y inside dé dn, the potential must be regarded as the 
sum of contributions from all smaller elements into 
which dé dy may be subdivided. 

Now the potential at a point P will be determined by 
conditions only at those points whose after-cones in- 
clude P—i.e., points lying in the Mach cone ahead of P, 
asin Fig. 2. We will consider a point P slightly above 
the element dé dn, so that the fore-cone of P intersects 
the element dé dn in a parabola, whose coordinates 
and £, are defined by 








do = (5) 


(6) 








B?[(y — m)? + 27] = (x — &)? (7) 
2 
y 
ra eS eps an 
Fic. 2. 


The element dé dyn will be further subdivided into 
smaller elements d£’ dyn’, in all of which the source in- 
tensity is constant at +g. The potential at P is then 
the integral over the area within the parabola of the 
contributions from these separate elements, 


2 Bs ™m 
¢ = -a f dé’ f 
a —-m 


ENE, .. AS 

V(x — #)? — Bly — 9’)? + 2°] 
This becomes 
AG — 2) |* 


t’)? cat B72?! =n: 





tf a V(x - 


Using Eq. (7) this is 


GOS YS leat Oi, Pe 
¢ = 2 f x dt’ = 3 Bz — a) 


and so 
0¢/0z = w= ng, or g= w/r (8) 


This result will hold even as z approaches zero and a 
approaches x, that is, for all points on the element dé dy. 
Since from Eq. (6) there is no contribution to the ver- 
tical velocity at P from points removed from this par- 
ticular element dé dn, it follows that the total vertical 
velocity over a source element defined by Eq. (5) is 
w= mg. The horizontal perturbation velocities at any 
point are affected by all source elements in the field 
and must be computed by an integration over the entire 
fore-cone of that point. 

Apparently, then, the boundary conditions for a thin 
airfoil will be satisfied by a distribution of source ele- 
ments of the type of Eq. (5), where w is a function of 
£ and » defined by the local slope of the airfoil according 
to Eq. (2). The pressure at any point on the airfoil is 
computed from Bernoulli’s equation, which, in this 
linearized form, becomes* 


P-fArA= —puU 


Here /p, is the free stream static pressure. A pressure 
coefficient may be defined as 


CE abt Pee (9) 
p 1/o9,U? U 





INFINITE WING, OBLIQUE TO FLOW 


Suppose a two-dimensional, infinite wing to be swept 
back at angle o to the flow. Two cases may be con- 
sidered, (1) tan o < 8, i.e., the leading edge lies ahead 
of the Mach wave from any‘point on it, or (2) tana > 8 
The first case is shown in Fig. 3. 


Bray ope 
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Fic. 3. 
‘ 

The pressure coefficient at a point P, (x, y), which for 
convenience may be taken on the x-axis, will be deter- 
mined by the source distribution within the fore-cone, 
A,of P. This forecone is bounded by the Mach lines 


rg n = +(x — £)/B (10) 


* See, for instance, Glauert (reference 3) 
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and by the leading edge 
£ = kn (11) 


where k = tane. The potential at P will be 


1 w dé dn _ 
ri + § j Ve-p or 


A new variable may be introduced to replace &, 


s=&—ky (13) 








It is seen that s is constant along lines parallel to the 
leading edge, and therefore the slope of the airfoil is a 
function only of s; 


\ = (w/U) = w(s)/U (14) 


Eq. (12) then becomes’ 


‘ 1 f 4 f- dn 
_. W as — 
® JO m W(x —s — kn)? — Bm? 


The limits of integration are, from Eq. (10), 
nm = —(x — £)/B, m= (x — §)/B 
or, using Eq. (13), 








: eo | 


oie 








2 = > n=- 


aos) 
+ 
>~ 


This can be rewritten as 
1 = m 
(= - = f w ds 4 
us 0 nz 
dn 


oa nme - — 

B(x — s)? | (x—s)k fam ey] 

a ee ee pees 5s ee ae ——————— — k2 
B? — k? Vee" 








This may be integrated to 


a 1 ( Lee ho) 
Qg = x Jo V/-k E B(x ei 5) | 


(2 


Or 
= w/(s) ds 
to- Jace 
The x-velocity increment is 
w(x) UX, 


=¢, = — a = SE 
wii V p? — k? Vp? — k? . 


The slope in the x-direction, \,, is related to the air- 
foil slope measured normal to the leading edge, A, by 


r 
A = Acoso = = 
V1 + ik? 
and, similarly, 
kr 


= —Asing = — = 
»y V1i+ k 


The pressure coefficient on a surface element will be 
C 2r, (16) 
nd 16 

P V 8? a 


This pressure produces a horizontal force on the air- 
foil in a direction normal to the leading edge, since the 
slope of the airfoil is everywhere normal to the leading 
edge. The magnitude of the horizontal force is given by 

2dr 2? 
= aes 8 SSS SCC”) 
Vet— ke Vite VB RP 

It will now be convenient to refer this drag coefficient 
to a dynamic pressure based on the velocity component 
normal to the leading edge, 


Cpu? 


Cy’ = = Cp(1 + k?) (18) 


u* cos? ¢ 
Furthermore, for the Mach Number normal to the 
leading edge, M, = M cos o we see that 
M? i B? — k? 


M,? -—1= a 
1 + k? 1 + k* 





If we define a new 8, referred to this normal Mach 
Number, as 


s’ = VM,2 —1 (19) 
we find that (for one surface of the airfoil) 
Cy’ = 2r7/p' (20) 


This is identical to the well-known result deduced 
by Ackeret for a two-dimensional supersonic airfoil, 
referred here to the velocity and Mach Number of the 
flow normal to the leading edge. The initial require- 
ment that 8 > k ensures that M, is greater than one. 
The effect of the sweepback is thus that the forces 
produced are determined only by the component of the 
flow normal to the leading edge—a result which can be 
seen quickly by direct physical reasoning and which 
was pointed out in 1935 by Busemann.® For a finite 
airfoil the result will, of course, be altered. 

A similar computation may also be made for the sec- 
ond infinite wing case, in which tana > 8. ‘This case 
is illustrated in Fig. 4. The-fore-cone of P will now 
intersect both leading and trailing edges of the wing, 
so the limits of integration must be revised somewhat 
The integration may be carried out in two parts, over 
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the areas bounded by the upper and lower Mach lines, 
respectively, from P. Using the new variable s we 


have again, 


Lf w ds dn 
QD = _ r = — a — 
T A V (x — s)? — 2(x — s)kn + (Rk? — B?)n? 


where A is the fore-cone of P. Writing the integral in 
two parts and transforming, 


¢=,+/, 


one obtains 


Be... a 
(x — sk we-s) g(x — s)? (21) 
\ Ee 5 ie lan k* — 8? 
h=— if w ds 
w os a 
dn 
(x i s)k , B2(s — x)? (22) 
5 a 
vive r ie alt Pi 


In the limits of integration, 7; and m2 are the coordi- 
nates of the Mach lines, 


Trak. ms s—k 
is k 4 ey 2 — k th 8 
and the lower limit, a, will be allowed to approach minus 
infinity. The first integral will be 


I, = f diaod cosh~—! k Ls | ae a4 x 
Io VR? — B? B B(x — s) Jle 


At the limit » = m, the expression in brackets becomes 


unity, so that 


1 2 (R? ech 2 
i= a d w ds cosh {2 — —- Pp 4 
rV k? — 8? Jo B B(x — t) 


The second integral is similar except that, in order to 
obtain the principal value, it must be noted that (s — x) 
is positive in the range considered. The expression in 
brackets under the square root is always positive in both 


cases. Therefore, 


| c (k? — B*)n)]|" 
nV k? — B? Js 8B p(s — x) J, 


or 


l “¢ k (k? — B")a 
I= — wdscosh-!| — — — 
rv k? — Bs B Bis —x) 


These integrals may be combined into 
l ra F ia k = (k® — B*)a 
$=- — w ds cosh—' |— — —— 
av k? — Bp S/o . B B(s — x) 
(23) 
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The incremental x-velocity is then 


Op eee 
Ox T 0 


aw(s) ds 
a 24) 


(x — s)V [a(x —s) — (k? — 6%)a]? — B(x — 5)? 


If a is now allowed to approach minus infinity, we 


Op _ I [ w(s) ds = 
Ox ed 2a 8B? 0 (x — s) =)) 


The total drag on the wing will be the integral over 
its surface of the product of pressure times local 


slope; the drag per unit span will be 


D= — > f A(x) Uo, dx 
0 


Using Eqs. (2) and (25) this is 


find 


ihe 
Ce 


p ‘ , w(s) ds 
D=-— — w(x) dx (26) 
aV k? — B? Jo 0 (x —s) 


This integral and _also the integral for @ in Eq. (23) 
are improper due to the singularity atx = s. How- 
ever, both can be shown to converge; moreover, the 
integral in Eq. (26) is identically zero, which can be 


shown by setting. 


p ¢ *wis) ds 
p= -— ————> w(x) dx + 
arV k? — B? 0 o (x —s) 
r “w(s) ds 
w(x) dx (27) 
0 x a = oy 


The second integral is over the triangular region 
shown in Fig. 5. Interchanging the order of integra- 


S 











Y 
x< 





Fic. 5. 


tion, and then interchanging variables, this second inte- 


gral becomes 


. “Sqw(x) dx " "* w(s) ds 
w(s) ds —— i — w(x) dx — 
0 0 (x— Ss) 0 0 (x—S) 


Thus, the second integral in Eq. (27) cancels the first, 
and the drag, as given by Eq. (26), is zero. 
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This result again can be anticipated from elementary 
considerations. In the present case, with k > 8, the 
Mach Number of the flow normal to the wing leading 
edge is less than one; this is equivalent to a two-di- 
mensional wing in a subsonic flow, with a superposed 
velocity parallel to the leading edge. Since this super- 
posed velocity should not affect the forces on the wing, 
the drag, as in any subsonic, perfect-fluid flow, should 
be zero. It is interesting to note, however, that in the 
case considered, an infinite body in a supersonic flow 
may have zero drag. For a finite wing, this will not 
be true. 


DRAG OF TRIANGULAR OR ‘‘DELTA”™’ WINGS 


It is desired to compute the form drag of an airfoil 
with plan form as shown in Fig. 6. A simple profile 
will be considered, consisting of the double wedge seen 
also in Fig. 6. The slope of the airfoil surface is con- 
stant in regions A and B, and the boundary condition 
can therefore be satisfied by a distribution of constant 
strength sources in A, and another in B. 


A ' 
TAN” K, 








The plan form is characterized by the sweepback 
angles of the leading edge, tan! R,, and of the maxi- 
mum thickness line, tan~' ko. In computation of the 
pressure distribution over the surface, three cases must 
be considered : 

(1) The Mach wave lies behind both the leading 
edge and the maximum thickness line, i.e., 8 > k; > 


(2) The Mach wave lies in between, i.e., kj > 8 > 


(3) The Mach wave lies ahead of both lines, 1.e., 
ki > ko > 8B. 

The fundamental problem involved is the computa- 
tion of the pressure distribution over a delta in which 
there is a constant source distribution. The finai 
wing may be obtained by superposing two such deltas. 
Two cases may then be distinguished: k > 8, andk < 
B. 





-—_> X& 


B(2 -y) 








Case (a): Single Delta,k >8 * 


In Fig. 7 it is desired to compute the pressure at point 
(x, y) due to a constant source distribution in the region 
bounded by the two Mach waves forward. from (x, y), 
and by the airfoil leading edges. The velocity potential 
will be the sum of three integrals 


od 0 «— Bly — 9) 
o(x,v) = — Lf nf F(&, 0) dE + 
v v1 kin 
y 2— B(y — 9) 
[ dn [ f(&, 0) dé + 
0 kin 
ye x + B(y — n) 
[ im [ f(é,n) dé (28) 
JY kin 


f(&, n) = [(x — &)? — By — n)?]7'” 


With proper attention to principal values, this be- 


where 


comes 


0 . 
w x + kn 
g=-— cosk—!—-— dn + 
7 FA B(y—n)" 
'y2 Paereee 
f cosh—' —~— ~ dy (29) 
0 Bin — y} 


Since the integrands vanish at y = y, and yo, the x- 
component of the disturbance velocity is 


- es dn 
ee + 
LJ, W(x + kn)? — By — 0)? 


‘v2 dn 
7 ee = (30) 
0 V(x — kn)? — B2(n — y)? 


This integral becomes 


a hoe : E as | + 


== Vk? — B ) B(x + ky) 


kx —. B*y |) 
OS a eee £ ‘ - 
Cc sh E bli aA ( $1 ) 


Since ¢, must be constant along lines x/y = constant, 
we introduce the variable ¢ = ky/x and also let n = 
k/B8. Itis seen that ¢ = 1 on the leading edge. Then, 
after some simplification, 
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2 ‘> (t/n)* (32) 


Ww 
=—-——=; cosh~' 
Oe BV nt 1 1-f 


This applies to0 <¢ <1,” > 1. 

In composition of the complete wing, the source 
distribution in B, (Fig. 6), may iufluence the pressures 
in some parts of A, so that the ¢, for points lying ahead 
of the delta in the present case must also be computed. 
After an analysis similar to that above, it is found that 








w 


2 ft. 9 
= — ——— cosh 33 
? ‘wera J — 


This applies tol <t<n,n > 1. 





Case (b): Single Delta, k < 8 


In this case, as seen in Fig. 8, the Mach wave from 
the nose of the delta lies behind the leading edge. In 
region A’, ahead of the Mach wave, the pressure must 
be constant and the same as that computed for the 
two-dimensional oblique airfoil, as given in Eq. (15), 
which may be rewritten as 


¢: = — Vina pager (34) 


+ wach wave 








Fic. 8. 


At a point (x, y) in A behind the Mach wave, the 
¢, may be considered to consist of the quantity given by 
Eq. (34), minus a contribution due to a source distri- 
bution in area C, Fig. 8. Evaluation of this contribu- 
tion results in 


ee a ae Ley) 
or Weal! 7 sin~" a (35) 


This is valid for0 <t<n,n <1. It may be noticed 
that ¢ = » on the Mach wave. 

These expressions for ¢, give essentially the pressure 
coefficient over a single-wedge wing with slope A, = 
w/U, since, from Eq. (9), 


C, = —2¢,/U 


Dp 


Hereafter, all slopes, A, will be measured in the x-z 
plane, so that the subscript x may be omitted. The 
results are plotted in Fig. 9, with C,8/) as a function of 
tand m. It may be noted that for m > 1, the pressure 
distribution is more nearly characteristic of a subsonic 
flow; this might be expected, since the component of 
the Mach Number normal to the leading edge is, in this 
case, subsonic. For m < 1, the constant pressure in 














» 


o of — J * 2 “6 
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the region ahead of the Mach wave becomes higher as 
n approaches 1, but it may be shown that the integral 
of the pressure over the entire delta for m < 1 is con- 
stant, independent of n. 


Using the basic Eqs. (32) to (35), the complete drag 
on a two-slope, or double-wedge, wing, such as that in 
Fig. 6, may be computed. If the slope in area A is \i, 
the source strength in that area must be w; = AU. 
If the slope in area B is \», the net source strength in B 
must be w = A2xU. However, it is convenient to com- 
pose the airfoil of superposed source distributions, each 
constant inside a complete delta back to the straight 
trailing edge. In this case 


Ww = UU and W2 = (Az —_ hu) U (36) 


CasE 1: B>k, > ke 


In the first case mentioned previously, 8 > k; > 
ko, i.e., the Mach wave lies behind both leading edge 
and peak line. Tlie drag will be composed of three 
parts: (1) the effect of w, acting over A, Fig. 10; (2) 
the effect of w, acting over B; and (3) the effect of we 
acting over B. Since the pressure coefficients are given 
by Eqs. (34) and (35), the first part of the drag will be 
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J. E — 2sintm, [t= Wn?) ia} (37) 
A" 7 1l- 


where A’ and A”, are the portions of area A shown in 
Fig. 10, and where the drag is computed for both sides 
of a symmetric airfoil. The quantity g is the dynamic 
pressure, 1/29U*. It is convenient to introduce a new 
parameter, 

ke No 1 ke 
=—=-, wheren, = —,m=— : 

ky ny Bi B 
The maximum thickness point is then located at r 
chords from the trailing edge. 

The elementary area appearing in Eq. (37) is then 
given by 


rT 


dt 


a 


dA = A7(1 — r)? 
where A; is the total area, Ay = A + B = k,b?/2 and 
where ¢ = kyy/x. The first part of the drag can then 
be written as 


Dy 8),2A r(l — r)? {f[ E 
= = Re| = — 
q BrV 1 — n" 0 2 


— l— (t/m)? dt \ 
m,| Py 1a — ir)? (39) 


The function of ” and r involved here must recur in the 
computation of the other parts of the total drag, so 
that it may be given a special name, 


i ae, 
Ginn) = Re {TO "|e 


rae 1 — (t/n)? dt \ ' 
sin nf] (1 ate tr)? (40) 
8r? 


Pu = A7Gi(m, r) 
q Bn 
The second part of the drag, due to the sources of 
strength w,, acting over B in which the slope is ¢, 
will be the difference of integrals taken over the entire 
area, Ay, and over the area A. Using the definition 
given in Eq. (40), this will be E 











Then 





(41) 


Dao SrAA2 


= —— A,[Gi(m, 0) — G,(n,, r)] (42) 
Br 

since G;(m, 0) will represent the integral of the pressures 

over the area A + B. 

The third contribution to the drag, due to the source 
distribution w,, can be represented by a similar func- 
tion, 

q Br 


BG, (m2, 0) (43) ° 


The expression G(n, r) can be evaluated by an inte- 
gration by parts and is found to be 





l—~? r 
Gi(n, r) = ore w= cos! ” + 


1 T aid 
en ae E + sin (rn) | (44) 


From this it is seen that Gi(m, 0) = 2/2. 

The two slopes of the airfoil must be functions only 
of the thickness of the airfoil and of the ratior. If +r 
is the total thickness expressed in chords (see Fig. 6), 
then the following identities may be easily proved: 


(7/2) [1/(1 — r)] 
— (17/2) (1/r) 
Mi — As = (1/2) [1/r(1 — 7)] 


Also B = rAr. Using Eqs. (44) and (45), the total 
drag coefficient, as given by the sum of Eqs. (41), (42), 
and (43), will be 
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where the subscription on m has been dropped. This 
expression is valid form < 1,and0 <r< 1. In the 
limiting case when m approaches zero, the drag coef- 
ficient becomes 


p = (1°/B) [1/r(1 — r)] (47) 


This is the correct result for a two-dimensional double- 
wedge airfoil, with peak located at r chords from the 
trailing (or leading) edge. Curves of Cp8/r? accord- 
ing to Eq. (46) are plotted in Fig. 11. 


CasE 2: 


In this case, the leading edge lies behind the Mach 
wave but the maximum thickness line does not (see 
Fig. 11). Therefore, none of the area A is influenced 
by the source distribution w, in B. The drag will be 
composed of the same three parts indicated in Case 1, 
although the contributions must be computed differ- 
ently. The pressure distribution over A will now be 
determined by Eq. (32). The drag contribution from 
this area due to the source distribution w, will be 


ki >B> kh 
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= => ~ ff cosh = WY a4 (48) 
@  prVn?—1 - # 


Using dA as given by Eq. (38), this becomes 
D 8A? ' 
—_ = a A, rG2(m, r) (49) 
q Br 


where 


G2(n, r) = 


(1 —r)? ; 
7 acai [. —————= COS —! 9 1 (t/ we a (50) 
Vn? —1 er 1l-f 


It can be seen that this is, in fact, the real part of the 
expression in Eq. (40). The second part of the drag, 
due to the influence of the sources w; in Ag, is 
Da Srdreo . . - 
asthe alas 71Ge(m, 0) — Go(m, r) | (51) 
q Bm 

The third contribution to the drag is again 

Pu, Bille — 

q Bar 
The integral in Eq. (50) can be evaluated, using 
an integration by parts. The expression for G: be- 
comes 





BG,(m, 0) (52) 
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V1 — rn? (— + LV? -—-l1- ‘y 


This expression is useful only for rn < 1, i.e., for r < 
1/n, which is the same as ky < 8. In the special case 
r = 0, it can be shown that 


G2(n, 0) = a. + sin-! : 54 
72(N, == a Ss jos (3 
‘ Vn? — 1 n sae 


The total drag coefficient, composed of the sum of 
Eqs. (49), (51), and (52), may now be written in the 
form 
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where the subscript on m has been dropped. This 


equation 1s useful only for m > 1, and rn < 1. 


CasE 3: ki > ke > B 


Both the leading edge and the maximum thickness 
lines now lie behind the Mach wave (see Fig. 13). We 
have now 2; > landrn, = m. > 1. The drag will now 
be composed of the three parts listed under Case 1, 
and an additional part due to the influence in A of the 
sources w: in B. This influence will be felt over the 
fraction A, of A, as shown in Fig. 13. 

The first contribution, D,, will again be computed as 
in Eq. (48) or Eq. (49), with the exception that the in- 
tegral defining G in Eq. (50) must be evaluated differ- 
ently, since 7m > 1. This might be denoted by 

Pu ae A7G2'(m, 7) (96) 
q Br 


where the prime denotes evaluation of Eq. (53) so as 
to be real for rn > 1. This may be written 


log n r cosh~‘ n 
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The second and third parts of the drag are again given 
by 
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SUPERSONIC WAVE 
Dao SAjA2 : oe? co 
= A[Ge(m, 0) 7 Ge (m1, r)] (98) 
g Bar 
and 
De SAo(Ao +. r - 
—= —— ) (99) 


BG2(nz, 0) 
q Br 

The new contribution to the drag due to the influence 

of we in A, will be the integral of the pressure defined by 

Eq. (33), in which the variable ¢ is now constant along 

rays through point M in Fig. 13, taken over area Aj. 


That is, 


Dy S02 — Ms) £ 2 — 
= a) cosh-'/"2— 144 (60) 
q BrV m2 —1JA, fp] 


In the area A, it can be shown that 


(1 — r)? 
dA = 8B _ = dt 
(f — r)* 
The drag contribution in A; due to w. can, therefore, 
be written 
Dy = 8A(A2 — Ad) : 
= — — BF (no, r) 


= (62) 
q Br 


where 


l—r)? sf 1 “peu 
/ ———— cosh-1.g/™—" ade (68) 
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The function F may be computed using an integration 


by parts, which results in 


F(n, 7) = 


a ¥ log n 
.F(n, r) = 4 | ae. Sd 
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For computational purposes it is convenient to replace 
n, by rm, and to define a new function by 
l1—rf logrn 


F’(n, r) = < 7 ——— 
l+rlvrn? — 1 


] vee? J (y2y? — 9 row 
ae log LV wr ~ tie = »|I (65) 


V n? n(1 — r) 
so that 
F'(m, r) = F(ne, 7) 
and 
Pa, Be = 3 ge. y) (66) 
q Br 


By inserting the proper values for \; and 2, the com- 
plete drag coefficient may then be expressed as the 
sum of Eqs. (56), (58), (59), and (66), This may be 


written 
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where the subscript on » has been dropped. 





(67) 


This expression is useful for m > 1, and rn > 1. 
However, the real part of Eq. (67) is the same as Eq. 
(55), so that the drag coefficient may be given by a 
single function. 

Using Eqs. (55) and (67), a family of curves of drag 
coefficient vs. r may be plotted with the sweepback 
ratio, n, as the parameter for values of m > 1. This 
is done in Fig. 14, using the abscissa (1 — 7), since the 
maximum thickness point is located (1 — r) chords 
from the leading edge. 

From inspection of Figs. 11 and 14 it is seen that the 
drag coefficient of the delta wing becomes higher than 
that of a two-dimensional wing when m < lL, i.e., as 
long as the leading edge is in front of the Mach wave, 
reaching a maximum when the leading edge is parallel 
to the Mach wave. Further sweepback will decrease 
the drag, but only appreciably if the line of maximum 
thickness is also behind the Mach wave. The drag 
will be less than the two-dimensional only when the 
tangent of the sweepback angle is greater than approxi- 
mately 1.4 times the cotangent of the Mach angle, 1.e., 


when 
k=tano > 14VM?—1 


Beyond this angle, the drag drops rapidly and is roughly 
half the two-dimensional drag when tan ¢@ is twice 














484 JOURNAL OF THE AERONAUTICAL SCIENCES—SEPTEMBER, 1946 


V M* — 1, if the maximum thickness -point is located 
at the proper place. 

It can be shown that as m approaches infinity, the 
drag coefficient approaches zero, at a rate given by 


os A log » (68) 
7? mr? on 

Therefore, if 8 is finite, i.e., the Mach Number is dif- 
ferent from 1, the drag goes to zero. However, if k 
is held finite, and the Mach Number approaches 1, 
then 8 becomes zero, and Cp becomes logarithmically 
infinite. However, the accuracy of the theory breaks 
down, of course, as this point is approached. 


Lirt OF A DELTA WING FOR 1 < 8 


In the case that the leading edge of a delta wing lies 
ahead of the Mach wave, i.e., m < 1, a solution to the 
problem of a flat plate wing at an angle of attack may 
also be obtained using the foregoing methods. In the 
upper half space, (g > 0) we may assume the solution 
produced by a uniform positive source distribution 
over the airfoil, resulting in a uniform upwash over the 
upper surface of the airfoil. In the lower half space, 
(z < 0), we may assume the solution produced by a 
uniform negative source distribution, resulting in a uni- 
form upwash also over the lower surface of the airfoil. 
The complete solution thus satisfies the boundary con- 
dition for a flat plate at a uniform angle of attack. 

This will be a correct solution, since, with the airfoil 
leading edge ahead of the Mach wave, no disturbance 
can exist ahead of the airfoil, and the boundary condi- 
tion is satisfied completely. This method of arriving 
at a solution cannot be used in the case that the leading 
edge is behind the Mach wave (i.e., 7 > 1) since the 
solutions for the upper and lower half space would pro- 
duce a pressure discontinuity in the plane z = 0, be- 
tween the airfoil leading edge and the Mach wave. The 
analysis for the latter case has been carried out by H. J. 
Stewart in a paper as yet unpublished. 

In the present case, the pressure distribution is anti- 
symmetric about the plane z = 0, with distributions 





over the airfoil surface given by Eqs. (35) and (9). 
The lift on both surfaces will be twice the integral of 
the pressure over the airfoil surface. But this integral 
is essentially that given by Eq. (39) or Eq. (41), with 
r = 0. Using Eq. (44), the surprising result is found 
that the total lift coefficient at angle of attack a is 
given by 


4a 4a 


6 ee egy (69 
6 VM -1 


But this is exactly the result for a two-dimensional, 
unswept-back supersonic airfoil. Thus, the lift coef 
ficient on a delta wing is independent of sweepback, 
and retains the value of the two-dimensional wing 
as long as the leading edge is ahead of the Mach 
wave. 


CONCLUSIONS 


Using the linearized equation for compressible fluid 
motion, solutions to the flow around any thin symmet 
ric body at zero lift can be found. As an example, 
the drag of a simple triangular or ‘“‘delta’”’ wing with a 
double-wedge profile has been computed. It is found 
that the drag coefficient of such an airfoil may be re 
duced to values much below that of a two-dimensional 
airfoil, but that this reduction will depend upon proper 
sweepback angle and location of maximum thickness 


point. The drag coefficient will no longer be a func- 
tion only of VM? — 1, as in the two-dimensional 
case. 


REFERENCES 


1 Ackeret, J., Luftkrafte auf Fltigel, die mit grossere als Schall 
geschwindigkeit bewegt werden, Z.F.M., Vol. 16, pp. 72-74, 1925 

2 von Karman, Th., and Moore, N. B., Resistance of Slender 
Bodies Moving with Supersonic Velocities, Transactions A.S.M.E., 
1932. : 

’ Glauert, H., The Effect of Compressibility on the Lift of Air- 
foils, Proc. Royal Society, (A), Vol. 118, p. 113, 1927. 

4 Prandtl, L., Theory of the Lifting Wing in a Compressible 
Medium, Luftfahrtforschung, Vol. 13, p. 313, 1936. 

5 Busemann, A., Aerodynamischer Auftrieb bei Uberschallge 
schwindigkeit, Proc. Volta Congress, 1935. 





In 
varia 
whic 
resul 
A sti 
quiré 
maxi 
not 1 
desig 
furtl 
prov 
spee 
air p 
requ 
whic 
supp 
of ai 
a me 
duct 
weig 


I" 


the 

forn 
bee! 
Ice: 
resp 
aerc 
thru 
pell 
occ 
tane 


Q). 

of 
Tal 
ith 
ind 


19) 


al, 


of - 


ng 








Icing Problems and the Thermal Anti-Icing 
System | 


F. L. BOEKE* ano R. A. PASELK* 


North American Aviation, Inc. 


SYNOPSIS 


In the design of the thermal anti-icing system the number of 
variables to be considered usually requires a lengthy analysis 
which, if the effect of the variables is not fully understood, may 
result in errors in heat required and unnecessary system weight. 
A study of the effect of altitude and air speed on the heat re- 
quired, pressure available, and pressure required shows that the 
maximum heat requirement occurs at sea level at some air speed 
not necessarily the maximum for the airplane. This maximum 
design speed is influenced by the adiabatic heating effect. A 
further study shows that in designing the ducting system to 
provide the proper airflow quantities, sea-level and minimum air- 
speed conditions should generally be used. The design of the 
air passages or gap adjacent to the surface to be heated usually 
requires acut-and-try procedure. In this paper a chart is shown 
which readily evaluates the effect of varying the gap size and the 
supply air temperature on the surface temperature, weight flow 
of air required, and gap pressure drop. Additional charts provide 
a means of selecting the optimum combination of gap and supply 
duct system pressure drop which results in a system of minimum 
weight and heat requirement. 


INTRODUCTION 


CE FORMATION ON AIRCRAFT has long been a hazard 
to commercial and military aviation as well as to 
the private flyer. Besides affecting the safety and per- 
formance of an aircraft, severe icing conditions have 
been the cause of a large number of flight interruptions. 
Ice formation on the wings, tail surfaces, and fuselage is 
responsible for an increase in weight and the loss of 
aerodynamic efficiency. Severe vibrations and loss in 
thrust may result from ice accumulation on the pro- 
pellers. During heavy icing conditions in which ice 
occurs on the propellers and the general airplane simul- 
taneously, the airplane may be unable to climb or even 
to maintain altitude and thus may be forced to land 
short of its destination. 

The problem of determining adequate and practical 
means of combating ice formation has therefore brought 
deicing and anti-icing systems to the forefront in air- 
craft design. 

This paper will deal with the theory of ice formation 
and the design of the thermal anti-icing system for air- 
craft, particularly as applied to the wing and tail sur- 
faces. Much of the information and theory used herein 
has been generalized from a review of the literature 
from 1928 to the present time. Most of the credit for 
the investigation and use of the thermal system may be 
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attributed to the National Advisory Committee for 
Aeronautics and the Air Technical Service Commandt 
of the Army Air Forces. 

Considerable progress has been made in this field 
within the last year but unfortunately most of the infor- 
mation is of a restricted or confidential nature and can- 
not be presented in this paper. 

An outline of the major points to be discussed is as 
follows: 

(1) Formation and Occurrence of Ice.—This is a study 
of the atmospheric conditions during which icing occurs 
and the theory of ice formation. It is presented as a 
logical introduction to the problem of thermal anti- 
icing system design. 

(2) Removal and Prevention of Ice.—A historical pres- 
entation is given of the methods used in the removal or 
prevention of ice and the heat requiiements of the 
thermal anti-icing system. 

(3) Design of the Thermal Anti-Icing System.—Time- 
saving methods are presented for the determination of 
the design altitude, air speed; heated air supply, double- 
skin air passages, and the ducting system. 


FORMATION AND OCCURRENCE OF ICE 


Icing Conditions 


In order that ice may form on an aircraft, water must 
be present in the atmosphere. Water occurs in the at- 
mosphere in the form of water vapor and in the form of 
water droplets. Two meteorological conditions during 
which ice may occur are a temperature inversion, and a 
cloud, fog, or mist. 

Normally, the air temperature decreases with alti- 
tude. However, a condition may occur where the 
temperature at a given altitude is higher than at some 
lower altitude. This condition is known as temperature 
inversion. For example, rain may be precipitating 
from a cloud that is above 32°F. and during its fall 
penetrate a strata in which the air is below 32°F. An 
airplane flying through this colder air strata may ex- 
perience icing due to the supercooled rain striking the 
cold airplane surfaces. 

In general, icing occurs in a cloud, fog, or mist when 
the temperature is below freezing, although ice may 
also form at temperatures slightly above freezing. 
Water droplets, which exist in the atmosphere at below 
32°F., are in a supercooled state and freeze upon colli- 
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sion with a surface. Concentration and size of these 
droplets, as well as the airplane speed and surface 
geometry, determine-the rate of freezing. Icing cases 
have occurred at temperatures as high as 36°F. and 
water droplets are believed to exist in the atmosphere 
at temperatures as low as —50°F. It is probable that 
the water droplet concentration and diameter decrease 
as the temperatures decrease, thus suggesting that the 
more severe icing conditions are likely to occur near 
32°F. 

Actual observations during icing conditions have 
indicated that the majority of serious ice accumulations 
occur at temperatures from 20° to 32°F., and seldom 
below O°F. or above 32°F. 


The Formation of Ice 


A number of factors, such as the airplane velocity, 
airplane geometry, and the prevailing atmospheric 
conditions, determine the type, severity, and extent of 
ice formation. 

Consider the paths of water droplets approaching an 
airfoil at high velocity. The water droplets do not 
follow the airflow lines but are deflected by them as they 
approach the airfoil, the amount of deflection being 
determined by the weight or mass of the water droplet. 
Thus, because of their greater weight, the large drop- 
lets will not be deflected so easily as the smaller drop- 
lets; consequently, more of them will strike the airfoil. 

Fig. 1 shows a plot of the paths of water droplets 
approaching an airfoil surface. The droplets at the 
outer edge of the region limited by the projected area of 
the airfoil in the direction of flow will not strike the air- 
foil, whereas those in the center, noted by the shaded 
area, strike the airfoil and, under conditions conducive 
to ice formation, will freeze. 


AIRFOIL PROJECTED AREA 
(AIR FLOW LINES «+ 





FLIGHT PATH 
OF AIRFOIL 





L__pATH OF A WATER DROPLET 
AS DEFLECTED BY THE 
AIR FLOW 


Path of water droplets approaching an airfoil showing 
deflection of water droplets by airflow lines 


ZONE OF WATER DROPLETS 
WHICH STRIKE THE AIRFOIL 


A supercooled water droplet is highly unstable and 
when it collides with the airfoil its supercooling is in- 
stantly canceled. Upon striking the airfoil the outer 
shell of the droplet evaporates, drawing its heat of evapo- 
ration from the center portion and causing that part 
of the droplet to freeze. This process is aided by the 
supercooling of the water droplet, convective cooling 
from the passing air, as well as by conduction to the air- 
foil. The rate of freezing is determined by the speed of 
evaporation, the ambient air temperature, and the mass 
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The greater the mass the longer the 
This difference in 


of the droplet. 
time required to remove its heat. 
the rate of freezing as determined by the water droplet 
size and air temperature automatically determines the 
type of ice formation. 

The water droplet and air particles approach the air- 
foil simultaneously. If the droplet is small it will 
freeze rapidly at the leading edge of the airfoil and the 
air particles have no opportunity to escape, resulting in 
a mixture of air and ice. This explains the phenomena 
of rime ice formation, which is white, opaque, and por- 
ous, its whiteness and porosity being attributed to en- 
trapped air particles. It also explains why rime ice 
always forms at the leading edge of an airfoil. 

Larger dropets take longer to freeze and will, upon 
striking the airfoil, freeze only in part at the leading 
edge, the unfrozen portion of the droplet freezing as 
it runs back over the airfoil. Since a longer time is 
required for freezing, air particles have an opportunity 
to escape, and clear ice, called glaze ice, forms along 
the entire surface. It is hard, dense, usually transpar- 
ent, and has good adhesive qualities. When frozen toa 
solid object, the bond between the ice and the object is 
stronger than the bond between particles of the ice and, 
consequently, becomes extremely difficult to remove 
except by melting. 

It may then be concluded that the formation of rime 
and glaze ice is fundamentally the same process, the 
difference in the type of ice being determined by droplet 
size and air temperature. 

A combination of rime and glaze ice is usually formed 
since various types and degrees of icing conditions may 
be encountered throughout the flight of an aircraft. 

Frost is formed by the condensation of water vapor 
into needlelike crystals of ice and occurs when the air- 
plane surface temperature is lower than the dew-point 
temperature of the surrounding air. The frost forma- 
tion is somewhat crusty but, having little strength, is 
unable to resist the high velocity airflow, changes in 
pressure, and vibration occurring in flight, and there- 
fore does not represent a serious icing condition Frost 
formation on an aircraft on the ground, however, re- 
quires careful consideration since it has been known to 
increase take-off distance and take-off air speed. 

The items that, when unprotected, may most affect 
the safety or performance of an airplane in icing condi- 
tions are: (1) carburetor; (2) propeller; (3) wing and 
tail surfaces; (4) pitot static tubes; (5) windshield; 
(6) antenna mast, loop, and antenna wire; (7) air 
scoops; and (8) engine cowling. 

These items are listed approximately in the order of 
their relative importance and in the order in which they 
should be considered. 

Although carburetor icing is of first importance, it is 
considered to be a separate subject and will not be dis- 
cussed further in this paper. 
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ICING 


Icing of propellers is next in importance since loss of 
propulsive efficiency due to heavy ice deposits and 
severe vibration directly affect the performance and 
safety of the aircraft. Propeller deicing is also con- 
sidered to be a separate subject and will not be dis- 
cussed in this paper. 

Heavy ice deposits on the wing and tail surfaces will 
result in loss of lift, as well as increased drag. The loss 
of lift, along with the increase in weight that accom- 
panies ice formation, will reduce the ability of the air- 
craft to ¢elimb or maintain altitude and in extremely 
severe conditions may result in loss of the aircraft. 

Ice on pitot static tubes, windshields, and radio 
antenna is a definite hazard and is the cause of air-speed 
and altitude errors, inadequate visibility, and poor radio 
communication. 

Icing on the carburetor air scoop, air intake scoops 
for the cabin heating system and the anti-icing system, 
and the oil cooler air intake scoop, affects the safety and 
performance of an aircraft. Probably of first impor- 
tance in scoop protection is the air inlet, which supplies 
air for the anti-icing system. 


REMOVAL AND PREVENTION OF ICE 


Numerous methods have been proposed and used for 
the protection of aircraft in icing conditions. These 
methods fall into two distinct categories—namely: de- 
icing or ice removal, and anti-icing or ice prevention. 

Water-soluble compounds, such as glycerin, honey, 
and Karo syrup have been used to lower the freezing 
point of water a few degrees below 32°F. This method 
was successful but required large quantities of the 
compound in order to provide anti-icing for a sufficient 
length of time. 

Various lubricants, in the form of greases and oils, 
have been used in an attempt to reduce the adhesive 
qualities of ice with respect to the aircraft surfaces. 
This method did not prove successful because the water 
droplets would imbed themselves into the lubricant 
layer and freeze, thereby actually accelerating the start 
of ice formation. 

The first practical means of protection consisted of 
an inflatable rubber boot attached to the leading edge 
of the surface. This system is designed to break the 
ice, once it has formed, by alternately expanding and 
contracting the rubber boot. Although this type of de- 
icing system is being widely used wherever protection 
against ice formation is required, its effectiveness in all 
types of ice appears to be limited. With increases in 
aircraft speeds, more difficulty will be experienced with 
keeping the boots attached to the airfoil surfaces, and 
the drag caused by the externally mounted boots will 
be more detrimental to aircraft performance and oper- 
ating economy. 

Some of the problems presented by the boot deicing 
system appear to have found a solution with the intro- 
duction of the thermal anti-icing system. The function 
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of the thermal system is to increase the temperature of 
the protected surface by means of heat. This can be 
done electrically or by means of circulating heated air 
within the surface. The hot-air system appears to be 
the most efficient at the present time from the stand- 
point of weight and power required, particularly for 
wings and tail surfaces. However, electrical systems 
may be more practical for protecting propellers and 


‘smaller items, such as pitot tubes, scoops, and radio 


masts. 

Heated air is usually obtained from heat exchangers 
installed in the engine exhaust system or from combus- 
tion heaters operating off the aircraft fuel system. 
Electrical power may be obtained from the normal 
engine-driven generators if the power requirement is 
small, or from auxiliary generators installed for the 
specific purpose. For each thermal anti-icing system, 
the designer must evaluate the available heat sources 
from a weight, power, and practical design standpoint. 
The design may show the use of a combination of heat 
sources, including electrical, engine exhaust, and com- 
bustion heaters, to be desirable. 


Heat Requirements for Anti-Icing 


Before any attempt can be made to design a thermal 
anti-icing system a knowledge of the amount of heat 
required to successfully remove or prevent ice is neces- 
sary. Early wind-tunnel tests were conducted by the 
N.A.C.A. to determine the rate of heat loss to dry air 
from airfoil surfaces. Initial flight tests conducted by 
the N.A.C.A. indicated that thermal anti-icing systems 
required that a 70° to 100°F. skin temperature rise in 
dry air be provided in order to give adequate protection 
in most icing conditions. In icing conditions during the 
winter of 1942-1943 the Army Air Forces in conjunc- 
tion with the N.A.C.A. tested a number of military air- 
craft equipped with thermal anti-icing systems. The 
results obtained.from these tests indicated the superi- 
ority of the thermal system and also showed that suc- 
cessful ice prevention could be obtained during a major- 
ity of icing conditions by transferring 900 to 1,200 
B.t.u. per hour per sq.ft. of the protected surface area. 

During this test period, however, a few severe icing 
conditions were encountered in which complete ice re- 
moval was not accomplished. This indicated that the 
heat to be transferred from the wing leading edge to the 
15 per cent chord point should be increased to provide 
better protection. 

Flight tests of the North American XB-25E heated- 
wing airplane at the Ice Research Base during the 
winter of 1943-1944 showed that approximately 2,000 
to 2,500 B.t.u. per hour per sq.ft. was being transferred 
from the leading edge area. This airplane has been 
particularly successful as an anti-icing system design 
because it was able to combat the most severe icing 
conditions at the Ice Research Base that winter. The 
heat transferred was considered more than adequate at 


that time. However, in the winter of 1944-1945 the 
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XB-25E encountered an icing condition, believed to 
be the heaviest experienced at the Ice Research Base 
during an icing test, in which it was unable to prevent 
run-back freezing although the leading edge surfaces 
were clear of ice at all times. This indicates that 
atmospheric conditions can be.extremely severe. 

Since the earlier tests, methods have been developed 
for determining the rate of heat transfer from airfoils 
which give good results for dry air. Recently much 
work has been done to determine the rate of loss from an 
airfoil section to wet air for various degrees of wetness. 
The heat requirement for this case is based on the 
amount of heat necessary to completely evaporate all 
moisture striking the leadingedge. The amount of heat 
to do this, however, must be determined by the atmos- 
pheric condition encountered. It varies with the con- 
centration of the water droplets in the atmosphere and 
the water droplet diameter, as well as the airplane veloc- 
ity and the shape of the airfoil. Preliminary calcula- 
tions by this method indicate that heat requirements 
may continue to increase at air speeds above 300 
m.p.h., which differs from the previous opinion that 
heat requirements would decrease above 300 m.p.h. 

There must be some upper practical limit to the 
maximum heat requirement for anti-icing. It is proba- 
ble that this optimum limit will be established when 
more data on the frequency of various severities of icing 
conditions are obtained. At the present time meteoro- 
logical data show that icing conditions would occur 
which, if considered in the design, would increase by 
four or more times the heat requirements now consid- 
ered satisfactory (2,000-2,500 B.t.u. per hour per sq ft.). 
Obviously, this is impractical. If such conditions do 
occur, they must be considered as occurring infre- 
quently and as being conditions for which protection 
cannot be provided. 

No study of the thermal anti-icing system is complete 
without mentioning the problem of ground defrosting. 
Frost formation on the ground is important since its 
presence on the wings requires increased air speed for 
take-off, as well as increased take-off distance. Not 
much is known about the heat requirements for ground 
defrosting but tests conducted on the XB-25E airplane 
indicate that sufficient heat is transferred to remove 
frost on the ground. It is believed that a system that 
will supply adequate anti-icing will also be adequate for 
ground defrosting, provided the air inlet scoop is located 
in a position where sufficient pressure recovery is avail- 
able from the propeller slipstream or provided some 
other power source is available. 


DESIGN OF THE THERMAL ANTI-ICING SYSTEM 


The thermal anti-icing system may be described as 
any design that is used for the prevention of ice by 
means of heat. However, this paper will consider only 


that type of thermal anti-icing system developed within 
the last few years for military aircraft. 


This design 


-obtain good heat transfer to the outer surface. 


utilizes heated air which is circulated through a narrow 
passage adjacent to the airfoil surface. In subsequent 
parts of this paper this heated air passage will be called 
the gap. As shown in Fig. 2, heated air is supplied toa 
spanwise duct from which it is discharged chordwise 
into the gap. The air, in passing through the gap, 
loses some of its heat to the surface and is then dis- 
charged through the aft portion of the wing to the 
atmosphere through outlet orifices. The purpose of a 
small gap is to create high air velocities and thereby 
How- 
ever, the minimum gap size is limited by the pressure 
available to the system as well as by production toler- 
ances. Extremely small gaps may be easily obstructed 
by metal particles, dirt, and other foreign material, 
which may eventually result in the failure of the system. 
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Fic. 2. Typical section through wing showing airflow path 
from spanwise duct into gap, path of air through spars, and air 
discharge. 


In order to establish a complete analysis of a thermal 
anti-icing system, the following must be considered: 

(1) The maximum surface temperature rise required, 
based on known atmospheric conditions or experimental 
test data obtained during icing flights. 

(2) The altitude and air speed at which the heat re- 
quired is a maximum and at which the ducting system is 
most critical. Also, the variation of heat required with 
altitude and air speed in order to determine the range of 
control required. 

(3) The maximum permissible supply air tempera- 
ture. 

(4) The most practical gap size and arrangement. 

(5) The design of the spanwise duct. 

(6) The design of the supply duct system. 

(7) Selection of the heater or heat exchanger. 

(8) Method of control. 


Determining the Critical Altitude and Air Speed 


The procedure of analysis required to establish the 
critical air speed and altitude for the heat requirement 
and the design of the ducting system generally leads to 
the investigation of a number of airplane operating con- 
dit ons. The procedure requires considerable time and 
may result in errors if the effect of velocity and altitude 
on the heat required is not known. The problem is 
approached by establishing a procedure that deter- 
mines the air speeds and altitudes which are critical for 
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the selection of the heat required and the ducting sys- 
tem. This will not only simplify the preliminary design 
calculations but the method is believed sufficiently 
accurate for use in the final design of the system. 

Before this procedure can be established, a study ot 
the following problems is necessary. 

(1) What is the heat required to prevent ice forma- 
tion as a function of airplane velocity and altitude? 

(2) What is the weight flow of air inside the gap as a 
function of the heat required? 

(3) What is the system pressure loss as a function of 
the weight flow of air? 

(4) What is the available air pressure differential as a 
function of velocity and altitude? 


For the Heat Required.—The heat transferred through 
the surface is 
Q = hAAt (1) 
where 


Q heat transfer rate 

h = heat transfer coefficient 

A area of heated surface 

At = surface temperature rise above ambient air 


ll 


ll 


It will be assumed that one selected maximum skin 
to ambient air temperature rise will be adequate for all 
air speeds and altitudes for a given icing condition. 
This assumption has been used in the design of most 
anti-icing systems to date. 

Since the area to be heated is constant for a given air- 
plane, the heat required per degree temperature differ- 
ence is directly proportional to the heat transfer co- 
efficient or: 


Q/At~h (2) 
~ = symbol meaning ‘‘is proportional to” 


Convective heat transfer from airfoils has been 
approximately correlated by the use of dimensionless 
moduli as: 


Nu = f(Re)”"(Pr)” (3) 
Nu = Nusselts Number = hl/k 
Re = Reynolds Number = /Vp/y = IG/u. 
Pr = Prandtl Number = C,u/k 


n,m = exponents determined by flow conditions 


Since the Prandtl Number is practically constant for 
air for the range of temperatures encountered in icing 
conditions, then 


Nu = f(Re)” (4) 


By substituting the variables contained in Nu and Re, 
Eq. (4) may be written 


hl ~~ ( vey’ or hl ~ (=) (5) 
k PA k rf = 


then 


k r n k : n 
h~ IVp orl ~ G (6) 
ES ay 


k = thermal conductivity of air 

} = any significant dimension of the airfoil 
V = true velocity of air over airfoil 

p = density of air 


= absolute viscosity of air 
G = weight flow per unit area 


The thermal conductivity k and absolute viscosity yu 
are functions of temperature only and can be assumed as 
constant for the range of air temperatures encountered 


in icing flights. Thus, Eq. (6) may be simplified to: 


h ~ (lVp)"/1 (7) 


A body moving in a high velocity air stream experi- 
ences a temperature rise due to adiabatic heating. The 
equation for the temperature rise At, may be written: 


At, = (0.85)1.8(V/100)? = 1.53(V/100)? (8) 


where V = true air speed in m.p.h. The correction 
factor, 0.85, has been shown by test to be a satisfactory 
modification in Eq. (8). 

The Aft to be used in Eq. (1) is then the required 
temperature difference between the surface and the 
free stream air, less the rise caused by the frictional 


heating. This is expressed as follows: 
At = Aft, — Af, (9) 
where 
At = temperature difference to be used in Eq. (1) 
for determining heat required 
At, = temperature of surface minus free stream air 
temperature 
At, = temperature rise of surface due to frictional 
heating 


To determine the effect of altitude and air speed on 
heat requirements, Eqs. (1), (7), and (9) are combined 


Q = K(A/Il)(IVp)"(At, — At,) (10) 


where K is a constant for a given flow condition or 
exponent 1. 

The effect of air speed on heat requirements may be 
determined by plotting Eq. (10) for constant altitude. 
This has been done in Fig. 3 for sea level, At, = 100°F., 
and for a fixed airfoil. Note that the heat requirement 
reaches a maximum at some particular air speed for each 
value of n. 

Fig. 4 shows the velocities corresponding to the 
maximum values of heat required shown in Fig. 3 for 
various values of At,, skin to free stream air tempera- 
ture differences. Fig. 4 was obtained by differentiation 
of Eq. (10) for various values of n. 

The exponent m may be obtained from flow transition 
data for the airfoil section used. In Fig. 5 an equiva- 
lent » value representing the average flow conditions 
over the area to be heated and the corresponding K 
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Fic. 5. Effect of transition point and heated chord length on 
average flow exponent. 


values to be used in Eq. (10) were calculated for various 
values of x/l, where x is the point of transition from 
laminar to turbulent flow measured downstream from 
the leading edge of the airfoil and / is the heated surface 
length. Values of m and K were approximated to ob- 
tain values of hk to within 3 per cent for values of 1Vp 
from 3 to 300. 
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| = heated length, ft. 
V = velocity, ft. per sec. 
p = density, lbs. per ft.* 


The transition point is a function of the following 
parameters: 


(1) Pressure distribution over the airfoil, particu- 
larly the location of the minimum pressure point. 

(2) Reynolds Number. 

(3) Mach Number. 

(4) Surface irregularities. 

(5) Air-stream turbulence. 


Little information is available on the position of the 
transition point, particularly with respect to its effect 
on heat transfer from airfoil sections. It is believed, 
however, that if the transition point can be determined 
from aerodynamic data, a constant distance can be 
added to it in the direction of airflow along the airfoil 
with favorable results for heat transfer calculations. 
The purpose of adding a constant length to the transi- 
tion point is to compensate for a gradual rather than 
abrupt heat transfer coefficient gradient which is be- 
lieved to exist as a result of the change from laminar to 
turbulent flow over the airfoil. 

Flight-test data on the XB-25E airplane indicate that 
if approximately 4 in. were added to the theoretical 
minimum pressure point, an effective transition point 
for the purpose of obtaining heat transfer coefficients 
would be established which would give the same average 
flow conditions and exponent m as was indicated by the 
measured heat transfer. 

Fig. 6 shows a plot of the effective transition points 
obtained by the above method for a high-speed airfoil. 
Transition data will vary considerably for different air- 
foils and therefore each airfoil must be individually 
considered. 

The step that follows is to combine the data in Fig. 6 
with Fig. 3 as shown by the curves in Fig. 7. The 
maximum heat requirement, as well as the air speed at 
which it occurs, is then known and is indicated by point 
X. 

General equations for the combined effect of altitude 
and air speed are established by taking ratios 
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Fic. 6. Plot of transition point versus indicated air speed. A 
plot of transition data for a high-speed airfoil section. 
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Fic. 7. Speed of maximum heat requirement for fixed transi- 
tion combining the transition data shown on Fig. 6 with the heat 
requirement curves on Fig. 3. 


- (72) (= wn =) (11) 
Q2 Vep2/ \ Ato — Ata 
Subscripts 1 and 2 represent any two conditions for 


which a comparison is desired. 
At a constant altitude 


M. PH. 


pi = pr 
Then, rewriting Eq. (11) and substituting Eq. (8) 


Qi _ (Ys ~ Lave) (12) 





Qe V2/ L At, — 1.53( V2/100)? 
At a constant true air speed 
Vi = Va 
Al,, = Ala, 
Then 
Q:/Q2.= (p1/p2)" (13) 


At a constant indicated air speed 
piVi? = p2V2? 
or 
Vi/ V2 = (p2/pi)'”* (14) 
Therefore, substituting Eq. (14) in Eq. (11) and 
rewriting 


Q _ (“YT Ai, — maiden} (15) 
On . At, — 1.53(V2/100)? 





It will be noted from Eq. (15) that for values of p/p. 
> 1, Q:1/Q2 will always be greater than 1.00, thus sea 
level is critical for the heat required at constant indi- 
cated air speed. 

For the Ducting System.—The heat transferred to the 
surface from the air in the gap may be expressed by two 
equations. 


Q = WCiVin — tout) (16) 
Q =hA Css. = beein) (17) 






W = weight flow of air in gap 

C, = specific heat of air 

ti, = supply air temperature entering gap 
tout = air temperature leaving gap 


t = mean gap air temperature 


ave. 


turin = Skin surface temperature 


The average gap air temperature may be expressed 
for this comparison as the arithmetic mean of /;, and 
tout OF 


lave. = (lin + tout) /2 (18) 
Combining Eqs. (16), (17), and (18) and solving for 
the heat transferred 


Q ie tin an sin (19) 
(1/hA) + (1/2WC,) 


Eq. (6) may be rewritten for flow inside the gap as 


n~= (78) or kw (20) (20) 
DX D\uA, 


W/A, 

cross-sectional area of gap 

significant gap dimension 

exponent determined by flow conditions in the 
gap 

k/u" varies only slightly with temperature (within 3 

per cent for f,ye. = 150° to 300°F.) and the gap size and 

arrangement is fixed; therefore 


h~ Ww” (21) 








~ 


II 


TORO 
ll 


or 


hA = CW” (22) 
C = a constant for the gap size and the gap air 
temperature 
Rewriting Eq. (19) . 
W*' (tin — tetcin) 


23 
(1/C) + (1/2C,W* ~ *) ar 


Q= 





In order to show the effect of the denominator in Eq. 
(23) on the heat required, Table 1 has been compiled 
for n’ = 0.8 for turbulent flow in the gap. 





TaBie 1 
% Variation 
1 1 - for Constant 
Values of C~ —2C,W®?~-C 2C,W®? Gap 
Gap Size—Inches 1/16 1/4 Wis 1/4 1/16 1/4 1/6 /, 


W,lbs. perhour10 1.68 6.83 1.32 1.32 3.00 8.15) 0% 
W, lbs. per hour 100 1.68 6.83 0.83 0.83 2.71 7.66f °“ 


6% 





This table shows that a variation of not over 10 per 
cent will occur for changes in W of 1,000 per cent, for 
the range of gap sizes generally used in thermal anti- 
icing system designs. Changes in gap temperature 
will not affect these results more than 3 per cent from 
tave. = 150° to 300°F. Thus, the expression (1/C) + 
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(1/2CpW' ~ "’) may be assumed constant for this com- 
parison, and since f¢,, and fi, have been assumed con- 
stant 
0:/Q2 = (Wi/W2)" (24) 
The general Fanning equation for pressure losses in 
straight duct may be expressed as 
AP = (pV?/2g)(4fL/D) (25) 
This equation may be rewritten as 
AP = C(W?/p) 26) 
where C is a constant and can be evaluated for any 
particular ducting system. Wis used instead of W7-5, 
which is correct for straight duct for turbulent flow 


since the largest part of the system pressure loss is 
assumed to occur in elbows and fittings where the loss is 


proportional to W?. 
From Eq. (26) for any two conditions 


AP, / AP, = (W,/W2)*(pe2/p1) (27) 
For a constant altitude 
pi = pe 
AP,/ AP, = (W,/Ws2)? (28) 


For the Pressure Differential Available.—The pressure 
available or dynamic head may be expressed as 


q = (pV?/2g) (29) 


For constant indicated air speed the pressure avail- 
able is constant. Since the lift coefficient and angle of 
attack are constant tor a constant indicated air speed, 
inlet scoop and discharge orifice pressures and ram re- 
covery efficiencies should be constant. 


Then, for constant indicated air speed 


qi/gz = 1 (30) 
and for constant altitude 
qi/ge = (Vi/ V2)? (31) 


Changes in indicated air speed are usually accom- 
panied by variations in ram efficiency and discharge 
orifice pressures due to change in angle of attack, and 
therefore the pressure available may not vary as shown 
in Eq. (31). However, at this point Eq. (31) will be 
assumed to be the correct relation and the variations 
with angle of attack will be discussed further in the 
conclusions. 

In any ducting system the pressure available must 
be equal to or greater than the pressure differential re- 
quired for a given airflow in order to obtain that airflow, 


or 


g/AP 2 1 


where g = pressure available and AP = pressure 
differential required for a given airflow. 


In order that a comparison may be made between any 
two altitudes or any two air speeds, the ratio of g/ AP 
will be assumed as equal to one for sea level, and for 
minimum indicated air speed. Thus, if for any other 
altitude or air speed g/ AP is less than one, the altitude 
or air speed being compared will be critical for the de- 
sign of the ducting system. When g/ AP is greater than 
one, sea level or minimum indicated air speed will be 
more critical than the altitude or air speed being com 
pared. 

Summarizing the equations developed for any two 
altitude conditions for constant indicated air speed 


a _ a Al Mh - peel ated be 
— = {| — — dD) 
At, — 1.53(V2/100)? 





Qs» pe 
Q:/Q2 = (Wi/W2)”’ (24) 
AP,/ AP, = (W1/W2)*(p2/p1) (27) 
g/g: = 1 (30) 


Rewriting Eq. (27) and raising the above expression 
to the n’/2 power, 


(ay = (my = Qi (32) 
\ AP 2p2 W2 Qz a, 


Combining Egs. (32) and (15) and rewriting 


AP, = ey oly ini 
AP, po K-—1 


where K = Afy/1.53(V2/100)?. 
At sea level for a given ducting system 


qi/ AP, = 1 or qi = AP, 
For constant indicated air speed g: = gz. Therefore, 


AP, = qi = qe 
and 


AP;/ AP; = q2/ AP; (34) 
Combining Eqs. (33) and (34) 


Ge (ey [ey (35) 
AP, pe K-1!1 ; 


A study of Eq. (35) reveals that when p;/p2. > | 
gz/ AP: can only in one case be less than one, that is 
when n/n’ < 1. 

However, g2/ AP2 can also be greater than, or equal 
to, one when n/n’ < 1, depending upon the numerical 
value of K, pi/pz and n’. In our previous discussion 
we assumed g;/ AP, = 1 for sea level, and when q,/ AP; 
< 1 altitude would be critical for the design of the duct 
system. In order to determine the altitude that would 
be critical, Fig. 8 was drawn for n’ = 0.8, ga/ APs = 1, 
At, = 100°F., and for various values of V and n. 

For each value of m in this figure, g@/AP:; = 1; for 
the region to the right of the given m curve q2/ AP; > 1; 
and for the region to the left of the given # curve 
g2/AP, < 1. Therefore, if the intersection point for 
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Fic. 8. Determining critical altitude for the duct system. A 
plot of Eq. (38). 


any given maximum icing altitude and minimum indi- 
cated air speed falls to the right of a given m curve, sea 
level will be critical for the design of the ducting system; 
if it falls to the left of the given m curve, the assumed 
maximum icing altitude will be critical. 

In the usual case, n > 0.75 at the low air speeds be- 
cause of early transition from laminar to turbulent flow 
at high angles of attack, the minimum air speed is 
greater than 150 1.A.S., and the maximum icing altitude 
is not over 20,000 ft.; thus, sea level is usually critical 
for the design of the ducting system. The example 
plotted in Fig. 7 shows n at 0.75 up to an air speed of 220 
m.p.h.1.A.S. Assuming that the minimum air speed = 
180 m.p.h. I.A.S. and the maximum icing altitude is at 
20,000 ft., and plotting these values in Fig. 8, the point 
of intersection falls to the right of the curve for » = 
0.75. Therefore, g,/AP: > 1, showing that sea level 
would be the critical altitude for the ducting system. 

Summarizing the equations developed for any two 
air speeds at constant altitude 


Qi a (2) Al, — 1.53( cm (12) 





Qs V2) LAt, — 1.53(V2/100)? 
Q:/Q2 = (Wi/W2)” (24) 
AP,/ AP, = (Wi/W2)? (28) 
gi/q2 = (Vi/ V2)? (31) 


Raising Eq. (28) to the n’/2 power 
(AP,/ AP2)"/? = (Wi/W2)" = Q:/Q2 ~—-(36) 
Combining Eqs. (36) and (12) , 


j 7\2n/n' hr Re R3 27]2/n’ 
SP, _ (, ' | 1.53(Vi/100) | (37) 





AP, V2 a 1.53( V2/100)? 
At minimum air speed for a given ducting system, 
AP, = qi. 


From Eq. (31) 
Qi = qo( Vi/ V2)? 


therefore 
AP, = a = go( Vi/ V2)? 
and 
AP; _ (FY) (38) 
AP, AP2\ V2 


Substituting Eq. (38) in Eq. (37) and rewriting 


ge SO naedelaas |: — 1.53(V; ory ¥ te 
—_— = -— —__-—- —— - (39) 
AP, \Ve At, — 1.53(V2/100)? : 

The following observations of Eq. (39) are made. 
When JV; represents a minimum air speed 


| — 1.53(V1 aor > 
At, — 1.53(V2/100)? 








If n/n’ S 1 then g2/ AP: > 1; if n/n’ > 1 then q2/ AP; 


may be less than 1. 
Eq. (39) is evaluated in Table 2 for values of n/n’ = 


(0.8/0.5) = 1.6. 








TABLE 2 
, Air Speeds M.P.H. = V;——— —_ 
Vi/V2 100 150 200 250 300 
0.8 0.8 0.83 0.888 0.976 1.12 
0.6 0.58 0.662 0.84 1.12 1.95 
0.4 0.47 0.76 Lae 7.67 
0.2 0.377 : af 





/ ° j a 
q:/APz2 vs. air speed (n/n'’ = 1.6). 


Note that many of the values of g2/ AP, < 1, particu 
laily at the lower air speeds, indicating that when the 
flow in the gap is partially laminar and the flow over the 
airfoil turbulent, speeds above minimum air speed may 
be critical for the ducting system. However, the flow in 
the gap is usually considered as fully turbulent. This 
effect will not be studied further but it does point out 
the fact that minimum air speed may not always be the 
air speed most critical for the duct system. 

A further point that should be considered is the 
effect of air speed on ram recovery efficiencies at the 
inlet scoop and the static pressures at the outlet orifice. 

In cases where g:/AP2 may appear to be sufficient 
for the higher air speed, a decrease in the ram recovery 
efficiency with increase in air speed may result in the 
pressure available being less than the pressure required; 
thus, the higher air speed may then be most critical for 
the design of the duct system. 

The conclusions to the study of the effect of air speed 
and altitude on the heater size and ducting system are 
summarized as follows: 

(1) The maximum heat requirement occurs at sea 
level for a constant indicated air speed. 

(2) At constant altitude, the air speed at which the 
heat requirement is a maximum may be determined 
with sufficient accuracy by the method given in this 
paper and as shown in Figs. 3, 4, 5, 6, and 7. 

(3) For most cases, sea level should be used for the 
design of the duct system because the pressure available 
may not be adequate for the weight flow and heat re- 
quirement at sea level. As a check, however, a graph 
for a constant required skin temperature to free stream 
air temperature rise as shown in Fig. 8 may be used to 
determine the critical altitude condition. 
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(4) For a constant altitude, minimum air speed may 
in most cases be considered as critical for the duct 
system. However, when the flow in the gap is par- 
tially laminar and the flow over the airfoil turbulent, or 
when ram recovery efficiencies decrease with increases 
in air speed, the higher air speeds may be more critical. 


Origination of the Design Chart for the Gap 


The engineer who has designed or aided in the design 
and calculations of thermal anti-icing systems is aware 
that it is difficult to analyze the heat balance in the gap 
without resorting to trial-and-error solutions or to 
approximations that may lead to the accumulation of 
errors. Further, the effects of gap sizes, supply tem- 
peratures, and gap pressure losses on heat requirements, 
as well as the balancing of flows between the upper and 
lower gaps, are sometimes difficult to visualize. There- 
fore, the authors have attempted to correlate the vari- 
ables used in the gap analysis on a single design chart as 
shown in Fig. 9. As the design criteria change a new 
chart can be made for the new requirements. 
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Fic. 9. Design chart for determination of gap size, gap pressure 
drop, and weight flow for a constant supply air temperature. 


Derivation of the actual formulas used in the con- 
struction of the chart is a lengthy procedure, conse- 
quently no attempt will be made to show their deriva- 
tion or the method of constructing the chart. How- 
ever, the time required to construct a chart of this type 
does not detract from its usefulness since one chart may 
be sufficiently general to apply to most designs. 

In the analysis of the gap the following equations 
must be solved 


Qi-s = WC, (ti — ts) (40) 
Qr-2 = he-3A (t, — ts) (41) 
Qs-4 = hs-4A(tzs — t4) (42) 
AP = (2fl/D)(pV?2/2g) (43) 


By proper combination of the above equations it was 
found that the weight flow, gap size, and pressure drop 
could each be expressed directly in terms of two para- 


meters. One parameter is a function of conditions out- 
side the airfoil, such as altitude, velocity, degree of 
turbulence, and heated area. The second parameter is 
a function of the air and surface temperatures. 

The equations for the parameters to be used in Fig. 9 
are as follows: 

For the gap, 


bl 
$9.8), 0.8 





K, = 0.131 (44) 


For the pressure drop, 


= 9g (bo 8b?-4\ / Py om 
Kap = 1.427 X 10 ‘( a \?) (45) 





For the weight flow of air in the gap, W 
Kw = ® = h3-4A (46) 


Q:-s = Qo-3 = Q;-4 heat transfer to the surface 
ho-3 = heat transfer coefficient, inside of gap 


h3s—4 = heat transfer coefficient, over outside of 
surface 

P = barometric pressure at altitude 

Po = barometric pressure at sea level 

W = total weight flow of air in the gap, upper or 
lower 

A = areaof heated surface 

V = velocity in the gap 

AP = pressure drop through gap 

l = length of gap, chordwise 

b = spanwise width of external heat transfer 
area = 1 ft. 


b, = spanwise width of gap, adjacent to 1 ft. of 
heat transfer area, ft. ; 


D = equivalent diameter of gap = 2d approx. 
d = gap thickness, ft. 

f | = flow friction factor in the gap (turbulent) 
K, Ka,, Kw = temperature parameters 

t; = supply air temperature 

te = log mean (f; — és) 

tj = surface temperature 

t, = ambient air temperature 

ts = air temperature leaving the gap 


Fig. 10 illustrates the location of the regions denoted 
by subscripts in the above equations as 1, 2, 3, 4, and 5, 
as well as the construction of the gap. 

The basic advantages for using this chart are: 

(1) The gap size and arrangement can be selected 
directly for a given AP and f3. 

(2) When two different surfaces are to be balanced 
for the same pressure losses, the gaps for both surfaces 
can automatically be selected. 

(3) If both the gap size and pressure loss are known, 
the skin temperature f; can be determined directly. 

(4) The weight flow of air can be read simultaneously 
with the other variables. 

(5) The chart provides a quick check on tests. 
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Lam al 
SECTION A-A 
DETAIL OF GAP CONSTRUCTION 


a DESIGN TEMPERATURES 


Fic. 10. Construction of gap and design temperatures showing 
location of region of heat transfer and significant gap dimensions. 


The design chart assumes that little or no heat trans- 
fer from the air in the spanwise duct to the skin surface 
occurs, which would normally increase the overall heat 
transfer efficiency. However, the supply air tempera- 
ture used in the chart can be varied to account for a 
spanwise temperature drop in the duct due to this heat 
transfer. 

Fig. 9 was used for the following studies of the effect 
of gap size and supply air temperature on gap pressure 
loss and the required weight flow of air. 

Effect of Gap Size on AP and W.—For a given flight 
condition and surface temperature, /;, the variation in 
the weight of air needed for heat transfer and the pres- 
sure loss may be determined for several gap sizes as 
shown in Fig. 13. 

Note the increase in pressure losses for small gap sizes 
and the reduction in weight flow of 75 per cent for a 
change in gap size from 0.250 to 0.040 in. 

Effect of t; on AP and W.—The curves of Fig. 11 show 
that increasing the supply air temperature results in a 
reduction of both weight flow and gap pressure drop. 
Approximately 35 per cent less heat is required for 
400°F. air than for 300°F. air. 

From this curve it is evident that a higher supply air 
temperature will provide a more efficient design. The 
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Fic. 11. Pressure loss and weight flow versus supply air tem- 
perature showing the effect of increasing supply air temperatures 
on decreasing heat requirements. 


maximum supply air temperature will be limited, how- 
ever, by structural overheating problems. 


The Optimum Combination of Gap Size and Duct System 


In order to determine the gap supply air duct and 
heater size, the designer usually allots a percentage of 
the pressure differential available to each of these items. 
This method of proportioning is based on judgment and 
past experience and will result in an adequate system 
from a pressure drop standpoint but it may not provide 
the minimum system weight and maximum system 
economy. 

The analysis that follows is unique in that the gap 
size and system duct diameters are determined simul- 
taneously, in order to obtain an optimum overall sys- 
tem. The method presumes that these data are avail- 
able to the engineer: 

(1) A general layout showing the airplane arrange- 
ment, so that the system equipment and ducting may be 
located and approximate length of ducting may be 
determined. 

(2) Critical air speeds and altitude as determined in 
an earlier part of this paper. 

(3) Surface areas to be heated and surface geometry. 

(4) Aerodynamic information on the airfoil section 
used, such as minimum pressure points, or other data 
that will establish airflow characteristics. 

(5) General equations for pressure loss in the ducting 
system. 

(6) General equations for heat transfer and pressure 
loss in the gap or a design chart similar to Fig. 9 in which 
these equations are evaluated. 

(7) Data showing weight flow vs. pressure drop for 
heaters of various types and capacities. 

A sample system, consisting of a combustion heater 
and ducting to wing and tail surfaces, is shown in Fig. 
12. Figs. 13, 14, and 15 present a solution for a system 
consisting of a heater and duct supplying air to the 
wings only. Figs. 16, 17, 18, and 19 present a solution 
for a system having a heater and duct branches to both 
the wing and tail surfaces. To illustrate the method a 


DUCT LENGTH INCLUDING 
EQUIVALENT LENGTH FOR 
HEATER =20 FT + 4 ELBOWS 
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Fic. 12. General layout for a duct system in an airplane. 
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Fic. 13. Gap size versus pressure loss and weight flow for the 
wing, for a particular system design, air speed, and altitude. 
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Fic. 14. Pressure loss versus weight flow and duct diameter for 
wing anti-icing system. 
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Fic. 15. 





Gap size and system duct diameter versus pressure 


loss for the wing system, showing points of minimum pressure 


loss for each duct diameter or gap size. 
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sample problem is worked out and is shown on the 
figures in heavy dotted lines. 

The critical air speed for this airplane is the minimum 
indicated air speed, and the critical altitude is sea level. 
The selected maximum supply air temperature is 350°F. 

The procedure for determining optimum gap size and 
ducting system is as follows: 

(1) The pressure loss in the gap and the weight flow 
of air for the required heat transfer vs. gap size are 
shown in Fig. 13. The data were obtained from the 
design chart (Fig. 9) and the heat requirement was 
determined for a given air speed, altitude, and surface 
geometry. In the case of a tapered wing, an average 
chord section should be used to determine an average 
weight flow of air, and the root chord section should be 
used to determine the gap pressure loss. 

(2) A plot of the duct system pressure losses fo1 
various duct diameters is shown in Fig. 14. The pres- 
sure losses may be calculated by means of accepted 
equations and expressed in terms of an equivalent length 
of pipe or number of diameters of pipe in the system. 
This may be determined from a general layout of the 
duct system as shown in Fig. 12. The weight flow of 
air to be used is the total weight flow of air required in 
the entire wing section. 

(3) The total pressure losses for the system are then 
determined by adding the pressure losses of the gap 
and the duct system for the same weight flows, and are 
plotted in Fig. 15. Note that for each duct diameter 
the pressure loss is a minimum at a particular gap size, 
which is the optimum for a given pressure available 
since the duct system is the lightest in weight. 

(4) For a given pressure available, a gap size and 
duct system combination may be determined for the tail 
surfaces. Figs. 16, 17, and 18 were obtained in the same 
manner as for the wing system. The selection of an 
optimum system was in this case limited by a minimum 
gap size determined by structural tolerances. 

(5) The duct systems for the wing and tail surfaces 
were each determined for the total length of duct from 
the air inlet to the surface. At point A in Fig. 12 the 
two duct systems join in a common duct. The equiva- 
lent length of this common duct may be determined and 
curves for several duct diameters drawn as in Fig. 19. 

(6) When the weight flows and duct diameters deter- 
mined for the wing and tail duct systems are plotted on 
Fig. 19 a horizontal line fixes the pressure drop in the 
common duct. By adding the two weight flows and 
plotting them on Fig. 19 at the same pressure drop, the 
common duct diameter is established. 

In the foregoing procedure the heater pressure loss 
was not considered. This may be included in the analy- 
sis by assuming an equivalent length of straight duct 
for the heater, or by subtracting the heater pressure loss 
from the total pressure differential available as used in 
the analysis. Both methods require that an estimate 
of the heater capacity and pressure drop be made. The 
latter procedure is preferred since a comparison of the 
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Fic. 19. Weight flow and duct diameter versus pressure loss; 
method of determining common duct diameter as a function of 
the wing and tail duct diameters and weight flows. 


effect of various types and sizes of heaters on the system 
can be obtained. 


CONCLUSIONS 


The following conclusions may be drawn from the 
studies of the gap: 

(1) A graphical procedure can be set up to simplify 
the gap analysis as well as to enable the engineer to 
visualize the effects of changes in the gap design on 
total system design. 

(2) The highest supply air temperature that can be 
safely used without overheating the structure should be 
considered, since it reduces the system weight and heat 
requirements. 

(3) A gap size can be selected which will result in an 
optimum system design, without resorting to tedious 
trial-and-error computations. The method shown 
may be enlarged upon to include more variables than 
are noted in this analysis; however, the method as 
presented will aid the designer in establishing optimum 
preliminary designs without excessive work. An im- 
portant advantage of the method is that all of the basic 
parts of the system may be established at one time, 
such as: the gap size and arrangement, the ducting 
sizes for all branch systems, and the heater size. 
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Chart for the Design of Wooden Spars 
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INTRODUCTION 


QO’ THE MANY TASKS with which an aircraft designer 
is daily faced, one of the most tiresome is the trial- 
and-error selection of structural sections, frequently 
requiring many repetitions of time-consuming compu- 
tations. In many cases it takes a considerable amount 
of labor to find a section which is just satisfactory, and 
a good real more to find the most efficient one. One 
such problem is the design of built-up spars for wooden 
aircraft. To assist in this case, a chart is presented by 
means of which allowable bending moments can be 
rapidly determined for spar sections of the types com- 
monly used. To aid in shear strength calculations, 
values of Q/J are also given. A section “‘efficiency”’ is 
defined on the basis of resistance to bending, and some 
typical graphs of this quantity are given. 


ALLOWABLE STRESSES 


Allowable stresses are based on the form factor for 
compressive stress developed by Newlin and Trayer' 
and are given in graphic form in reference 2, Fig. 2-6. 
The only limitation on the computed tensile stress sug- 
gested in reference 1 is that it shall not exceed two and 
one-half times the maximum computed compressive 
stress. Many designers, however, prefer to set an 
arbitrary limit on the tensile stress, somewhat lower than 
what is implied by the recommendation of reference 1. 
In the preparation of the chart given here, this limit 
has been set at the modulus of rupture of solid rectangu- 
lar sections, which for spruce is 9,400 Ibs. per sq.in. 
This point is not of great practical importance, however, 
since in most cases the inverse load conditions on the 
wing will ensure adequate lower flanges. 


NOTATION 


F, = compression modulus of rupture, lbs. per sq.in. 

I = moment of inertia of section about neutral axis, in.‘ 
Q = first moment of area about neutral axis, cu.in. 

4 = areaof section (booms plus webs), sq.in. 

S, = I/(h — y), section modulus, compression side, cu.in. 


5S = I/y, section modulus, tension side, cu.in. 

M, = allowable bending moment, compression critical, in.Ibs. 
M, = allowable bending moment, tension critical, in-Ibs. 

G = asection parameter 

y = 5'/b 


The meaning of other symbols is illustrated in Figs. 
and 2. 


—" 
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FORMULAS USED 


The area assumed to be effective in bending is that 
shown hatched in Fig. 1. For normal thin web spars, 
this approximation introduces but small error. It 
gives nearly the same result as assuming, when there 
are two webs, that only one is effective in bending, 
which is fairly common practice. (When the modulus 
of elasticity of the webs differs from that of the booms, 
the designer may wish to use an “‘effective’’ web thick- 
ness in determining 0.) Using this assumption, the 
reader will easily verify that the section properties can 
be expressed in nondimensional form as functions of 
t./h and k by the following formulas: 


a/h=[1/(1+k)]—'/2(t./h) 


a 1 ~ Saw 
h 1+k 2h 


h GP [(a/h)(t./h)k]/{G[1 — (5/h))} 


S,/bh?=G 
S,/bh?=G{ [1—(5/h)]/(5/h)} 


where 


1 t.\3 t./a\? 
—(1+?){ k(1+k)-{ - 
12 * (*) waNtT, iG) (1) 


G= { iti —o i 
1- (¥, /h) 


For spruce spars, using a limiting tensile stress of 9,400 
Ibs. per sq.in., the allowabie bending moment is then 
given by the lesser of 


M,/bh?= F.G (2) 
or 


M,,/bh?=9,400G { [1 — (5/h)]/(9/h)} (3) 


a 
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ine) 


Fie. 2. 


When M,.=M,, i.e., when the allowable stresses in 
tension and compression are reached simultaneously, 


then 
F,=9,400{ [1 —(9/h)]/(9/k)} (4) 


DESCRIPTION OF CHART 


In the lower right quadrant, solutions of Eq. (1) are 
given. Values of k up to 2.4 have been included to 
cover the case of normal spars subjected to inverse 
loading. The compression modulus of rupture, F,, is 
given in the upper left quadrant. ‘This portion of the 
chart is‘reproduced from reference 2, Fig. 2-6. 

In the upper right quadrant, the solid lines give values 
of M/bh?, and the dotted lines represent Eq. (4) for 
different k. These give the upper limit of M/bh? for 





Design chart for spruce spars. y = b’/b. F, = 


compression modulus of rupture—lbs. per sq.in. + 1,000. 


the givenkandG. A value of F, greater than that cor- 
responding to the dotted line for the relevant k indicates 
that tension is critical. If, following reference 1, the 
designer does not wish to limit the tension stress in 
this way, the dotted lines should be ignored. 

Values of h(Q/J) are given in the lower left quadrant, 
which is used independently of the others. 


EXAMPLE 


A typical solution is indicated on the chart for a 
spar having t,/h = 0.22, y = 0.18, k = 0.6. Following 
the dashed lines in the directions of the arrows leads to 
the point A, with M/bh? = 900. Since A lies below 
the dotted line for k = 0.6, then compression is critical, 
and the maximum moment is given by M, = 900bdh’. 
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If the tension flange is reduced in depth to make 
k = 0.5, then the point corresponding to A falls above 
the line for Rk = 0.5. Tension is therefore critical, and 
the solution is given by the point B, with maximum 
moment M, = 800bdh?. 


OTHER MATERIALS 


When woods other than spruce are used, the allow- 
able moments obtained from the chart should be multi- 
plied by the ratio of the modulus of rupture in static 
bending of the species, as given in reference 2, Table 
2-1, to that of spruce. 

SECTION EFFICIENCY 


It is suggested that a criterion for section efficiency 
be the allowable moment divided by the area, i.e., 


n = M/A in.Ibs. per sq.in. 
The area may be written 
A/bh = ¥ + (t./h)(1 + k)(1 — ¥) 


whence, 
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Fic. 3. Typical spar efficiencies. 


To illustrate the effect of the three variables on 
n/h, graphs of this quantity are given in Fig. 3, which 
shows typical variations with y, k, andt,/h. The dotted 
line gives the values when the tensile stress is limited to 
9,400 Ibs. per sq.in. The corresponding solid portion is 
for unlimited tension. Where there is only one line, 
the two values are the same. 
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jected to Transverse Loading Only, N.A.C.A. T.R. No. 181, 
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2 ANC-18, Handbook on the Design of Wood Aircraft Struc- 
tures. 
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